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PREFACE 


This  document  covers  research  as  to  the  feasibility  of  multlspectral 
mosaic  Image  conversion  panels  carried  out  at  the  G.  T.  Schjeldahl  Company, 
Inc.  during  the  twelve  months  from  1  February  1971  to  >1  January  1972. 

The  research  was  directed  by  Dr.  Donald  E.  Anderson;  co-investigator 
was  Richard  L.  Swisher.  The  research  staff  during  all  or  part  of  the 
period  covered  by  this  report  included  Robert  Mracek,  Eugene  Hildreth, 

Lloyd  Schultz  and  John  Rooks. 


ABSTRACT 


A  Development  program  is  described  in  which  the  feasibility  of  mosaic 
EL/PC  Image  conversion  panels  sensitive  to  UV,  near  IR,  X-rays,  and 
visible  light  was  studied.  Thin  Film  photoconductors  are  electroded  in 
a  regular  array  with  unit  cells  0.020"  on  centers.  These  arrays 
are  connected  to  opaque  electrode  arrays  forming  the  back  pads  of  an 
electroluminescent  (EL)  lamp  array  through  the  use  of  microglass  spacer 
sheets.  Both  thick  film  EL  and  thin  film  (TFEL)  lamp  arrays  were  prepared 
and  studied.  Useful  near  IR  response  was  obtained  for  inputs  of  10~6 
w/cml.  Useful  visible  light  response  was  obtained  for  5  x  10"^  ft-cd 
input.  Peak  optical  gains  of  500  were  achieved  with  thick  film  EL  and 
optical  gains  of  3000  using  TFEL  output. 

The  combinations  of  materials  used  were  prepared  in  test  sample  form, 
electrically  and  optically  parameterized,  and  then  computer  simulations 
were  performed  to  determine  the  range  of  parameters  needed  for  a  success¬ 
ful  assembly.  The  computer  models  simulate  the  transient  or  steady  state 
optical  stimulation  of  EL/PC  cells  with  simple  sinusoidal  power  applied 
or  more  complicated  waveforms.  All  computer  programs  used  are  docui^nted . 
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1.0  General  Introduction 


1 . 1  General  Description  of  Construction  and  Operation  of  EL/PC 

Image  Conversion  Panels 

A  solid  state  electroluminescent-phococonductor  (EL/PC)  image 
converter  is  a  device  which  responds  to  radiation  of  a  certain  intensity 
and  wavelength  distribution  by  emitting  radiation  of  higher  intensity 
and/or  different  spectral  distribution.  Our  studies  have  concentrated  on 
image  conversion  panels  having  a  structural  geometry  consisting  of  layers 
of  mosaic  patterned  cells.  All  photoconductors  studied  were  deposited 
as  thin  films  with  mosaics  of  patterned  electrodes.  Input  radiation 
is  incident  normal  to  the  PC  film  and  photocurrent  flow  is  in  the  plane 
of  the  PC  film.  Figure  1  is  a  back-lighted  photomicrograph  of  a  mosaic 
electroded  PC  layer  fabricated  in  our  laboratory.  The  opaque  material 
is  the  molybdenum  electrodes . 

To  connect  the  small  square  electrodes  of  the  PC  layer  in 
Figure  1  to  individual  EL  lamps  in  a  layer  above  the  PC  layer,  we  need  an 
interface  layer.  Ore  type  of  successful  Interface  layer  that  we  have 
previously  developed  is  shown  in  Figure  2.  In  this  photomicrograph,  the 
round  dots  are  "posts"  of  solder  alloy.  Each  of  these  posts  are  bonded 
to  one  of  the  small  square  electrodes  of  the  PC  layer  in  Figure  1.  The 
light  colored  matrix  is  0.003"  thick  glass.  This  is  used  for  dimensional 
stability  of  the  sheet. 


The  forward  or  upper  ends  of  the  conductive  posts  are  connected 
to  an  array  of  rear  electrode  pads  for  the  EL  output  layer.  If  thick  film 
EL  is  used,  it  is  sprayed  onto  the  electroded  interface  layer  and  a  cont¬ 
inuous  semi-transparent  front  electrode  is  applied.  Figure  3  is  a  photo¬ 
micrograph  of  the  light  emitted  by  an  array  of  thick  film  EL  lamps  fabri¬ 
cated  in  this  manner.  An  opaque  dielectric  layer  is  usually  applied  behind 
the  output  EL  layer  to  prevent  light  feedback  to  the  PC  layer.  If  thin 
film  EL  (TFEL)  is  used  as  an  output  media,  then  it  must  be  prepared  on  a 
separate  substrate  and  electroded  and  then  bonded  to  the  interface  layer. 

Figure  4  is  a  vertically  exaggerated  drawing  of  a  section  of  an 
EL/PC  image  conversion  panel  using  thick  film  EL.  The  operation  can  be 
visualized  by  tracing  the  flow  of  electrical  power  from  the  conductive 
lattice  across  and  modulated  by  a  stimulated  PC  cell,  up  a  conductive 
pos t ,  and  by  means  of  capacitance,  through  the  EL  layer  to  the  continuous 
front  electrode.  Figure  5  is  a  schematic  of  how  one  EL/PC  cell  operates. 

As  radiation  impinges  upon  the  PC,  its  resistance  decreases.  This  switches 
an  increasing  fraction  of  the  applied  voltage  onto  the  EL  capacitor.  In 
response  to  this  voltage,  the  EL  material  emits  light.  The  image  conversion 
panel  can  be  thought  of  as  an  array  of  photoresistive  voltage  dividers. 
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Extremely  Simplified  Schematic  of  Operation 

of 

EL/PC  Image  Conversion  Panel 


FIGURE  5 


1.2  Philosophy  of  Thl9  Program 


Under  a  previous  contract  (N00017-70-C-0213) ,  it  was  shown 
that  mosaic  EL/PC  infrared  image  converters  could  be  simulated  by  computer 
modeling  and  the  operating  characteristics  of  a  completed  test  assembly 
agreed  with  the  predictions.  This  type  of  simulation  makes  possible  many 
economies  in  the  search  for  useful  combinations  of  photoconductors  and 
electroluminescent  materials. 

In  this  program,  we  were  to  continue  studies  of  infrared  (IR) 
photoconductors  used  with  thick  film  EL  and  to  extend  this  work  to  photo¬ 
conductors  sensitive  to  ultra-violet,  visible,  and  X-ray  radiation.  Also, 
we  were  to  further  investigate  the  possibility  of  using  thin  film  EL(TFEL) 
in  image  conversion  panels.  If  the  computer  simulations  indicated  that 
a  combination  of  materials  could  function  properly,  then  test  assemblies 
were  to  be  fabricated,  if  possible,  and  compared  with  the  simulations. 

In  addition,  we  were  to  investigate  the  possibility  of  placing  a  visible- 
to-visible  image  conversion  panel  on  top  of  any  other  type  of  image 
conversion  panel  to  boost  the  conversion  gain  of  the  pair. 

The  work  and  this  report  was  divided  into  five  tasks: 

.  Develop  Thin  Film  EL  lamp  arrays  compatible  with 
possible  assembly  techniques 

•  Investigate  UV  sensitive  photoconductors  to  find 
possible  candidates  for  use  in  test  assemblies 

.  Investigate  X-ray  sensitive  photoconductors  to  find 
possible  candidates  for  use  in  test  assemblies 

.  Investigate  PbS-TFEL  combinations  and  improve  the  CdSe 
device  fabricated  under  previous  contract. 

,  Improve  CdS,Se  visible-to-visible  image  conversion 
panels  and  investigate  possibility  of  cascading  with 
similar  panels  or  any  of  the  image  conversion  panels 
above . 

2.0  Thin  Film  EL  Fabrication  and  Testing 
2.1  Introduction 


Thin  film  EL  (TFEL)  has  these  advantages  over  thick  film  EL: 
.  Higher  brightness 

.  Better  brightness  maintenance  characteristics 

•  More  uniform  appearance.  It  is  microscopically  uniform 
in  brightness  within  a  cell. 


TFEL  has  these  major  disadvantages: 


.  It  is  much  more  difficult  to  fabricate  consistent  lamps 
from  sample  to  sample.  Many  vacuum  deposition  processes 
are  involved.  The  heat  treatments  necessary  to  form  a 
good  lamp  limit  the  possible  assembly  procedures  of  an 
EL/PC  device. 

.  Its  brightness  is  a  very  steep  function  of  applied 
voltage  which  results  in  very  little  gray  scale  range 
in  an  EL/PC  image  conversion  panel. 

.  Its  capacitance  per  unit  area  is  much  larger  than 
thick  film  EL  and  hence,  its  impedance  per  unit  area 
is  lower.  This  requires  that  PC  cells  be  fabricated 
with  a  lower  range  of  impedances  than  those  necessary 
for  thick  film  EL. 

Despite  the  above  shortcomings  and  difficulties,  TFEL  is  still  a 
desirable  light  output  media  because  of  its  high  brightness  and  excellent 
life  if  made  properly.  These  characteristics  are  especially  important 
for  aircraft  displays. 

During  this  program,  we  improved  upon  the  performance  obtained 
from  TFEL  fabricated  under  the  previous  contract. 

2.2  TFEL  Fabrication  Techniques  and  Results 

2.2.1  Co-Planar  Electrodes 


It  was  originally  planned  that  we  try  to  increase  the 
impedance  of  the  TFEL  cells  by  using  a  co-planar  electroding  scheme  as 
shown  in  Figure  6.  Experimental  trials  indicated  that  such  a  configuration 
held  little  promise  for  success  because  of  extreme  demands  on  electrode 
edge  smoothness  and  electrode  spacing  of  less  than  0.001".  This  electrode 
spacing  is  impractical  for  large  area  panels. 

2.2.2  Series  TFEL  Lamps 

Another  impedance  raising  approach  was  then  tried.  This 
electroding  technique  is  shown  in  Figure  7.  The  technique  used  a  "moated" 
rear  electrode  for  the  TFEL  cell  and  a  "floating"  semitransparent  front 
electrode.  If  the  area  of  the  region  common  to  the  power  grid  and  floating 
electrode  and  the  area  common  to  the  "island"  and  the  floating  electrode 
are  equal,  then  the  impedance  of  the  TFEL  lamp  is  quadrupled  over  a  simple 
sandwich  electrode  structure.  The  total  operating  voltage  of  the  TFEL  lamp 
is  doubled  but  the  electric  fields  in  the  material  remain  the  same  as  for 
a  simple  sandwich  cell.  We  abandoned  this  approach  because  the  lamp  could 
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CO-PLANAR  TFEL  ELECTRODES 
FIGURE  6 


"MOATED”  DRIVE  ELECTRODES 
FOR 

TFEL  OUTPUT  CELL 


FIGURE  7 
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not  be  fabricated  as  in  Figure  7.  Because  the  TFEL  needs  high  temperature 
baking,  the  front  semitransparent  electrode  has  to  be  Sn02  or  In02  deposited 
on  a  substrate  first.  The  moated  rear  electrode  had  to  be  deposited  last 
and  we  could  not  find  an  etching  technique  for  the  metals  used  (Au,  Ni, 
or  Al)  that  didn’t  damage  the  lamps  in  some  unknown  way.  All  lamp  samples 
suffered  from  voltage  breakdown  problems. 

At  this  time  we  decided  to  concentrate  on  the  simple 
sandwich  cell  of  TFEL  and  try  harder  to  adjust  the  impedance  range  of 
the  photoconductors. 

2.2.3  Sandwich  Style  TFEL  Cells 

Concurrent  with  the  unsuccessful  attempts  at  using  novel 
TFEL  cell  geometries,  we  were  attempting  to  improve  the  performance  of 
"simple"  sandwich  style  TFEL  cells.  Figure  8  is  a  cross-section  through 
such  a  cell.  The  transparent  substrates  we  used  were  typically  boro- 
silicate  glass.  We  tried  two  types  of  semitransparent  front  electrodes: 

Sn02  and  In02.  The  phosphor  was  ZnS  doped  with  manganese.  The  dielectric 
is  GeOj.  The  rear  electrode  is  typically  Au. 

At  the  beginning  of  the  program,  we  performed  an  abbre¬ 
viated  matrix  study  using: 

a)  Three  types  of  Sn02  coated  substrates. 

b)  Three  types  of  Ge02  evaporation  sources. 

c)  Two  types  of  ZnS:Mn  phosphor  evaporation  sources. 

The  plan  was  to  abbreviate  the  matrix  by  building  lamps 
using  all  types  of  processes  in  (a)  above  with  a  standard  set  of  processes 
from  (b)  and  (c)  above.  The  best  (a)  process  was  then  used  with  all  types 
of  (b)  processes  and  the  same  (c)  standard.  This  optimization  was  carried 
through  each  successive  layer  of  the  lamp  structure.  The  types  of  substrates 
Included  commercially  obtained  substrates  and  Sn02  and  In02  coated  substrates 
prepared  in  our  laboratory. 

All  ZnS:Mn  depositions  were  made  in  a  bell  jar  vacuum 
system  using  an  internal  chimney  deposition  system  as  shown  in  Figure  9. 

T'.'.e  best  lamp  obtained  by  this  matrix  of  permutations 
was  fabricated  of: 

1.  Commercially  obtained  Sn02  coated  glass  substrate. 

2.  250 Jl  of  GeO  2;  vacuum  baked  2  hrs.  at  500°C. 

,  18,000^  of  ZnS  (shot  from  charge  containing  207. 

MnS);  vacuum  baked  2  hrs.  at  500°c. 
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FIGURE  8 
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The  AI2O3  Crucible  is  Heated  by  a  Tungsten  Coil. 

A  Close  Fitting  Plug  of  thick  Quartz  Wool  is 

Placed  Over  the  Deposition  Charge  and  Held  in 
Place  by  a  Simple  Refractory  Metal  Split  Ring. 


Diagram  of  ZnS,  CdSe,  pbS,  PbSe 
'  Deposition  Fixture 


FIGURE  9 
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4.  9,000A  of  GeOa. 

TM 

5.  Karma  rear  electrode. 

This  lamp  would  emit  1,000  ft-L  when  driven  at  260  V  rms  and  5KHz.  Life¬ 
times  of  this  type  of  lamp  were  only  ~50  hrs.  at  which  time  voltage 
breakdown  at  small  spots  would  begin  to  destroy  the  test  cells.  Decay 
was  very  slight  before  failure.  The  voltage  breakdown  problem  plagued 
us  throughout  the  program. 

Figure  10  is  a  representative  sample  of  the  brightness 
vs.  voltage  characteristics  of  the  TFEL  lamps  fabricated  in  this  matrix 
study.  The  results  are  typically  that: 

1.  Increasing  thickness  of  Ge02  layer  shifts  the  curves 
to  increased  operating  voltage. 

2.  Increasing  Mn  concentration  in  the  phosphor  layer 
(No. 3  above)  shifts  the  curve  to  lower  operating 
voltage  and  steepens  the  slope  of  the  curve. 

3.  Post  deposition  vacuum  baking  of  the  phosphor 
layer  is  very  important  for  maximum  phosphor 
brightness  and  lifetime. 

Figure  11  is  a  comparison  of  a  lamp  fabricated  in  the 
matrix  study  with  a  commercially  available  TFEL  lamp  (Sigmatron,  Inc., 
Goleta,  California).  Our  lamp  was  fabricated  as  follows: 

1.  Commercial  Sn02  coated  glass  substrate 

2.  250k  Ge02,  unbaked 

3.  15,000^  ZnS  phosphor,  shot  from  charge  containing 
20%  MnS,  vacuum  baked  2  hrs.  at  300  C. 

4.  9,000^  Ge02 

5.  Au  rear  electrode 

All  dielectrig  and  phosphor  depositions  were  made  with 
a  substrate  temperature  of  200°C,  The  ZnS  evaporation  crucible  has  a 
quartz  "wool"  plug  placed  on  top  of  the  charge.  A  thermocouple  encap¬ 
sulated  in  quartz  is  used  to  monitor  the  temperature  of  the  quartz  "wool" 
plug.  For  this  lamp,  the  temperature  was  1150  C.  The  lamp  operated  for 
240  hours  at  near  50  ft-L  before  voltage  breakdown  occurred. 

The  operating  frequency  was  5  KHz  and  the  voltage  was 
209  V  rms.  The  voltage  was  gradually  Increased  to  218  V  to  compensate 
for  aging.  All  good  lamps  exhibited  excellent  brightness  maintenance. 

But  at  some  time  after  approximately  cne  hundred  hours,  they  suffered  a 
catastrophic  breakdown  which  obliterated  the  entire  cell. 
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BRIGHTNESS,  FT-L 


Driving  Voltage,  rms  Volts  @  5 KHz 

BRIGHTNESS  VS.  VOLTAGE  OF  LABORATORY  SAMPLE  LAMP  COMPARED  WITH  AVAILABLE 

TFEL  LAMPS 

FIGURE  11 
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This  catastrophic  failure  mode  of  the  TFEL  lamps  which 
plagued  us  was  traced  to  almost  sub-microscopic  pinholes  in  the  phosphor 
layer.  Very  small  defects  in  the  surface  of  the  finished  lamps  were 
observed  to  be  the  center  of  destruction  of  macroscopic  areas  that  broke 
down  in  operation.  We  were  already  using  a  quartz  "wool"  plug  in  our 
phosphor  deposition  boat.  If  we  increased  the  thickness  of  this  plug 
or  used  a  less  porous  plug,  we  expected  deposition  times  to  become  very 
lengthy.  Consequently,  we  began  experimenting  with  thicker  nucleating 
layers  under  the  phosphor  layer.  Previously,  if  a  GeOs  nucleating  layer 
of  greater  than  ~250A  was  used,  the  lamp  surface  crazed  during  vacuum 
baking.  We  changed  the  Ge02  evaporation  source  from  an  A1203  crucible 
to  a  Pt  resistance  heated  source  of  simple  folded  ribbon  geometry.  Much 
better  adhesion  of  Ge02  to  the  front  electrode  was  then  observed. 

The  catastrophic  failure  mode  of  the  TFEL  lamps  was 
reduced.  The  deposition  source  for  the  Ge02  layers  and  the  increased 
thickness  of  the  phosphor  nucleating  layer  helped. 

The  best  lamp  fabricated  since  using  a  thicker  nucleat¬ 
ing  layer  was  manufactured  in  cross  grid  matrix  style  with  0.016"  wide 
In02  lines  on  0.020"  centers  for  the  semitransparent  electrode. 

The  lamp  was  fabricated  of: 

1.  In02  coated  substrate  (etched  into  lines). 

2.  500A  Ge02  unbaked,  deposited  from  Pt  boat.  Charge 
is  fused  before  opening  shutter. 

3.  20,000A  of  ZnS,  shot  from  a  charge  containing  207. 

MnS,  vacuum  baked  two  hours  at  500°C.  Source 
diffuser  plug  temperature  is  1150°C, 

4.  6,000^  of  Ge0z  deposited  from  Pt  boat.  Charge  is 
fused  before  opening  shutter. 

5*  Al  rear  electrode  (deposited  as  lines). 

This  lamp  was  operated  in  a  constant  brightness  mode. 

It  exhibited  a  fast  decay  mode  after  initial  turn-on.  Its  beginning 
operating  parameters  were  185  V  rms  at  5  KHz  to  emit  50  ft-L.  Its  "first" 
half  life  was  ~10  hrs.  At  this  time,  the  voltage  was  raised  to  192  V  rms. 
Its  projected  half  life  at  this  voltage  was  -  1125  hrs.  After  800  hrs. 
running  time,  the  lamp  required  200  V  to  emit  50  ft-L,  Small  portions 
of  the  lamp  did  suffer  catastrophic  failures;  final  failure  occurred  at 
approximately  1200  hours.  We  feel  that  this  particular  phosphor  fabrication 
technique  does  yield  a  useful  phosphor  layer  but  our  substrate  preparation 
procedures  need  further  refining.  At  this  stage  of  development,  the  TFEL 
effort  was  then  directed  toward  methods  of  integrating  TFEL  cells  into 
EL/PC  test  assemblies. 
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2.3  Interesting  Anomalies  of  TFEL 

While  fabricating  a  test  EL/FC  assembly,  a  small  portion  of  the 
assembly  was  coated  with  the  following  layers: 

1.  Glass  substrate 

2.  Molybdenum  base  electrode 

3.  ~  10,000l  of  CdS,  Se : C 1 

4.  3,OOoA  GeOa 

5.  16,000A  of  ZnS:Mn  (shot  from  a  charge  containing  15%  MnS) 

6.  3, 000 A  Ge02 

TM 

7.  Semitransparent  (~  40%)  Karma  electrode 

This  combination  of  materials  produced  a  very  low-voltage  TFEL 
lamp.  The  operating  characteristics  of  the  lamp  are  shown  in  Figure  12. 
The  lamp  would  not  emit  light  under  direct  current  stimulation.  As 
circumstances  permitted,  we  used  scrap  samples  of  TFEL  lamps  to  try  to 
determine  what  caused  this  behavior.  Our  best  guess  as  a  result  of 
several  loosely  controlled  experiments  is  that  Cl  doping  from  the  CdS, 

Se  layer  caused  the  phenomena.  As  a  check,  we  ran  a  ZnS.’Mn  TFEL  layer 
through  a  photoconductor  Cl  doping  cycle.  After  depositing  the  dielectric 
layer,  we  obtained  a  TFEL  lamp  having  operating  characteristics  as  shown 
in  Figure  13.  This  was  the  lowest  voltage  lamp  we  fabricated  (when 
operated  at  low  frequencies).  The  lifetime  of  the  lamp  was  only  a  few 
hours . 
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Light  Out,  Ft-L. 


Volts  rms 

TFEL  Lamp  with  CdS  Layer  Included 


FIGURE  12 


BRIGHTNESS 


3.0  UV-to-Vialble  Image  Conversion.  Panels 


3.1  Introduction 


We  originally  proposed  that  ultra  violet-to-visible  mosaic 
image  converters  might  be  feasible  if  large  band  gap  photoconductor 
materials,  those  greater  than  3  ev,  exist  which  are  compatible  with  current 
EL/PC  construction  techniques  and  exhibit  fast  photocurrent  rise  and 
decay  times.  Such  a  candidate  material  is  ZnO  (band  gap  ~  3.2  ev) .  The 
literature  shows  it  to  have  a  slow  photocurrent  rise  and  decay  time. 
However,  these  reports  do  not  clearly  indicate  its  behavior  in  thin  film 
form  under  the  large  signal  voltage  biases  necessary  for  mosaic  EL/PC 
image  conversion.  Thus,  materials  such  as  Al2S3,  ZnO,  ZhS,  ln203,  Sn02  , 
Sb203,  Bi203,  and/or  others,  as  indicated  by  brief  literature  search, 
wore  to  he  prepared  in  thin  film  form  and  doped  according  to  the  best 
information  available.  Photocurrent  under  near  UV  stimulation  was  to  be 
measured  while  ac  voltages  of  50  to  250  volts  rms  at  400  to  4,000  Hz 
were  applied.  These  measurements  of  photocurrent  when  used  in  the  computer 
simulation  of  EL/PC  devices  would  help  determine  the  feasibility  of  a 
UV-to-visible  image  converter. 

3.2  Results 


The  literature  search  turned  up  deposition  techniques  only  for 
ZnO  and  of  course,  ZnS.  ZnS  has  a  very  low  conductivity.  The  few  samples 
of  pure  ZnS  that  we  prepared  with  0.002"  electrode  spacings  and  thicknesses 
of  6,000A  to  20,000A  had  resistance  values  >  10^  Q.  we  could  not  detect 
any  photoconductivity  when  stimulated  with  intense  .336ti  radiation. 

We  then  turned  to  ZnO  as  it  is  somewhat  better  documented  as  a 
photoconductor.  It  seems  that  most  techniques  for  producing  ZnO  as  a 
photoconductor  are  held  as  proprietary  and  little  has  been  published  on 
the  subject.  We  first  attempted  to  deposit  ZnO  by  reactive  sputtering  in  a 
partial  atmosphere  of  02.  All  attempts  had  very  poor  adhesion  to  the 
substrates.  Next  we  used  r-f  sputtering  and  obtained  clear,  well  adhering 
films.  In  thicknesses  of  6,000A  to  18, 000 A  the  films  had  a  cell  resistance 
(.002"  electrode  spacing,  see  Figure  1)  of  1,5  x  105ft  to  lO^i).  However, 
we  could  detect  no  photoconductivity  under  intense  UV  stimulation.  We 
attempted  to  Cu  dope  the  films  by  imbedding  the  substrates  in  doped  powder 
and  baking  the  container  at  400  C  for  various  lengths  of  time.  The  only 
result  was  cloudy,  peeling  films.  Again,  no  photoconductivity  was  observed 
but  the  results  are  not  conclusive  since  the  films  were  very  crazed  and 
cracked. 


As  one  last  experiment  we  checked  to  see  if  the  ZnS:Mn  TFEL  lamps 
would  emit  light  at  less  than  their  normal  voltage  under  UV  stimulation. 

No  effect  other  than  simple  fluorescence  was  observed. 

As  a  consequence  of  the  above  negative  results  and  lacking 
information  as  to  proper  ways  to  proceed  in  post  deposition  treatment  of 
ZnO  films,  we  discontinued  effort  on  this  task. 
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4.0  X-Ray  to  Visible  Image  Conversion  Panels 


4.1  Introduction 


Past  vork  on  EL/PC  mosaic  structures  showed  that  uncooled  PbS 
photoconductlve  ceels  appeared  to  not  have  a  useful  dynamic  range  of  re¬ 
sistance  for  use  in  a  PC -EL  enfrared  image  converter.  Very  little  has 
been  done  to  investigate  its  use  as  an  X-ray  sensitive  photoconductor,  and 
it  is  judged  that  carrier  lifetime  is  probably  too  short  for  this  application. 

However,  a  hetero junction  device  of  PbS  and  CdS  or  PbSe  and 
CdSe  shown  in  Figure  14  was  proposed  as  a  possible  X-ray  sensitive 
photoconductor  combination. 

A  sound  theoretical  justification  for  expecting  useful  per¬ 
formance  from  such  a  device  couldn't  be  presented.  The  photoconductivity 
mechanisms  for  PbS  are  not  well  understood.  However,  work  performed  on 
a  proprietary  X-ray  dosimeter  utilizing  thin  films  of  Pb  as  X-ray  to 
electron  converters  indicated  that  such  a  heterojunction  device  should 
be  investigated. 

4.2  Results 


The  first  problem  encountered  with  the  structure  of  Figure  14 
is  that  we  could  not  prepare  PbS  with  small  enough  dark  conductivity  to 
adequately  turn  off  an  EL  cell. 

The  PbS  layer  would  effectively  short  out  the  CdS  layer.  The 
only  way  we  could  decrease  the  conductivity  of  the  PbS  layer  was  to 
decrease  the  thickness  to  ~  1000-2000A  and  sinter  it  in  an  02  atmosphere 
for  several  hours.  However,  films  that  were  this  thin  were  obviously 
too  thin  to  absorb  enough  X-ray  photons  to  be  useful.  The  only  other 
thing  we  could  try  was  to  see  what  happened  if  the  PbS  layer  was  patterned 
so  that  it  did  not  contact  the  electrodes.  If  a  high  resistance  layer 
formed  at  the  PbS-CdS  junction,  the  CdS  film  could  retain  a  smooth  enough 
electric  field  across  its  electrode  gap  to  function  as  a  photoconductor. 
Experimental  trials  indicated  that  the  opposite  effect  occurred.  It 
appeared  that  the  PbS-CdS  interface  became  a  very  conductive  region 
compared  to  either  material.  As  an  increasing  voltage  was  applied  to 
the  electrodes,  the  CdS  layer  would  begin  to  break  down  near  the  electrodes. 
It  acted  as  if  the  PbS-CdS  interface  was  a  current  shunt. 

At  this  point  in  the  program,  it  was  decided  to  cease  effort 
on  this  task  as  proposed. 

We  checked  to  see  If  the  operating  voltage  of  a  TFEL  lamp  was 
reduced  by  X-ray  stimulation.  There  was  no  observable  effect. 

4.3  Recommendations 


Using  thin  film  photoconductors  to  convert  X-ray  photons  into 
photocurrcnts  of  a  magnitude  suitable  for  controlling  the  emission  of 
electroluminescent  lamps  is  most  probably  not  practical. 
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HETEROJUNCTION  CONFIGURATION  FOR 
CARRIER  INJECTION  TRANSVERSE  TO  E  FIELD 


As  long  ago  as  eighteen  years,  moderately  successful  X-ray- 
to-visible  conversion  panels  were  fabricated  of  think  film  EL  and 
thick  film  CdS  layers.  These  devices  were  never  practical  due  to  construction 
problems.  Thick  film  CdS  is  very  difficult  to  prepare  in  uniform  layers 
of  the  necessary  ~  0.010'*  thickness  for  adequate  X**ray  absorption.  Also, 
if  the  devices  are  prepared  in  a  "sandwich"  structure,  which  was  the  only 
practical  process  in  the  past,  the  shunting  capacitance  across  the  PC 
limits  the  useful  range  of  the  device.  However,  if  TFEL  instead  of  thick 
film  EL  could  be  used,  the  effective  shunting  capacitance  across  the  CdS 
would  be  reduced  significantly.  If  further  work  with  direct  conversion 
of  X-rays  is  to  be  pursued,  the  TFEL  and  thick  film  CdS  or  CdSe  sandwich 
structure  should  be  investigated. 

It  should  be  noted  that  indirect  conversion,  using  a  conventional 
Patterson  CB-2  phosphor  screen  to  convert  X-rays  to  visible  light  and 
then  a  visible-to-visible  image  intensifier  panel  to  amplify  the  image, 
has  been  successfully  developed  under  another  contract*.  In  this  way, 
reductions  in  radiation  levels  have  been  achieved  which  are  quite  significant. 
Optical  gains  in  excess  of  100  were  achieved  on  10"  x  10"  panels  with 
0.020"  center-to-center  cell  spacing;  image  saturation  can  be  obtained 
with  less  than  20  mR  exposure  to  70  KeV  X-rays,  xhis  then  constitutes 
a  performance  "yardstick"  against  which  any  direct  conversion  approach 
must  be  measured. 


*U.S.  Army  Electronics  Command,  Fort  Monmouth, 
Contract  No.  DAAB07-71-C-0267. 
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5.0  Visible-to-Visible  Image  Conversion  Panels 


5.1  Introduction 


Under  this  part  of  the  program,  we  were  to  develop,  if  possible, 
two  basic  types  of  visible  light  amplifier  panels.  The  first  type  was 
to  be  a  very  sensitive  panel.  This  type  was  to  be  sensitive  to  as  small 
an  input  light  level  as  possible.  The  second  type  was  to  have  as  large 
an  output  brightness  as  possible  with  its  sensitivity  threshold  possibly 
no  better  than  10*2  to  10"!  ft-cd.  The  high  brightness  panel  was  to  be 
considered  mostly  for  use  as  the  second  stage  of  an  image  conversion  system. 
The  first  stage  of  such  a  system  could  be  whatever  type  of  image  converter 
we  developed  (IR,  UV,  X-ray,  or  Visible  Light  sensitive).  This  high 
brightness  second  stage  would  be  used  as  a  "brightener"  to  overcome  embient 
light  problems  in  viewing  an  image  on  the  first  stage  panel  by  itself. 

To  build  such  a  high  brightness  panel,  a  TFEL  output  is  virtually  required 
if  we  were  to  have  a  chance  of  achieving  many  tens  of  ft-Lamberts  output 
brightness. 


The  first  type  or  sensitive  panel  would  probably  have  to  be 
one  using  thick  film  EL  for  output.  Based  on  previous  experience,  thick 
film  EL  cells  are  a  much  better  impedance  match  for  sensitive  CdS  type 
photoconductors  than  are  TFEL  cells.  Also,  TFEL  cells  had  not  been 
successfully  integrated  into  an  assembly. 

As  was  mentioned  in  Section  4.3  above,  concurrent  with  this 
contract  we  were  also  working  on  a  contract  from  Fort  Monmouth  (UAAB07-71- 
C-0267)  to  develop  our  existing  light  amplifier  technology  into  the 
capability  for  fabricating  10"  x  10"  image  conversion  panels  using  0.020" 
EL/PC  cells.  This  excellent  research  mix  involving  an  assembly  technique 
development  program  greatly  benefited  both  programs.  We  used  the  PC 
doping  results  of  this  program  in  fabricating  10"  x  10"  light  amplifier 
panels  for  the  Army  and  we  used  the  improved  assembly  techniques  in 
fabricating  thick  film  EL  devices  in  this  program. 

5.2  Sensitive  Light  Amplifiers 
5.2.1  PC  Materials  Study 

We  began  with  three  types  of  PC  deposition  material. 

The  first  type  was  CdS.  In  all  doping  experiments  with  C'u  and  Cl  as 
dopants,  the  resulting  films  were  sensitive  to  light  levels  of  5  x  10*3 
ft-cd  or  greater  but  response  times  were  very  sluggish.  Response  time  for 
a  CdS:Cu,Cl  cell  to  reach  907.  of  its  final  conductivity  was  typically 
minutes  for  a  stimulation  of  5  x  10*3  ft-cd.  Also,  this  material  was 
very  history-sensitive.  After  an  hour  of  light  stimulation  for  testing, 
the  PC  cells  typically  had  to  be  kept  in  the  dark  overnight  to  be  able  to 
duplicate  the  results  of  the  testing. 


We  also  tried  using  a  mixture  of  80%  CdS  and  20%  CdSe 
as  a  deposition  material.  When  doped  with  Cu  and  Cl,  films  of  this 
material  exhibit  no  long-term  "history"  effects.  However,  it  was  difficult 
to  include  enough  Cu  in  the  films.  The  only  way  we  could  Cu  dope  with 
any  success  was  to  imbed  the  films  in  a  crucible  filled  with  an  overdoped 
(Cu2S)  CdS  powder  and  fire  the  covered  crucible  in  a  N2  atmosphere  at 
400°C  for  times  of  15  minutes  to  2  hrs .  This  method  would  increase  the 
resistivity  and  sensitivity  of  the  films.  But  the  films  are  cloudy  and 
have  poor  uniformity.  Subsequent  doping  with  HC1  gas  further  reduced  the 
uniformity  but  sensitivity  was  generally  increased. 


With  the  poor  results  of  the  previous  two  types  of  PC, 
we  returned  to  the  material  we  have  used  for  several  years.  This  is  an 
80%  CdS,  20%  CdSe  finely  divided  powder  already  doped  with  Cu  and  Cl  and 
quite  photoconductive.  Before  depositing  this  material,  it  must  be  baked 
in  vacuum  at  600°C  for  at  least  sixty  minutes  to  drive  out  the  chlorides. 
The  PC  films  will  not  adhere  to  the  substrates  if  this  is  not  done.  By 
using  this  doped  powder,  the  PC  film  received  adequate  Cu  doping  during 
formation.  Prior  to  this  contract,  we  had  attempted  to  Cu  dope  the  PC 
films  by  evaporating  minute  amounts  of  Cu  during  PC  film  deposition  but 
the  results  were  always  erratic. 

Figure  15  is  a  plot  of  the  PC  cell  conductivities  of  a 
2"  x  2"  PC  array  made  with  the  pre-doped  CdS,Se  powder.  The  bars 
represent  the  ranges  of  conductivity  observed  by  measuring  25  cells  spaced 
uniformly  about  the  substrate.  The  straight  line  drawn  through  the  plot 
is  the  least  squares  fit  to  the  data  used  in  the  simulation  programs  dis¬ 
cussed  in  the  Appendices. 


This  type  of  photoconductor  can  have  its  impedance 
values  at  a  stimulation  level  of  10“*  ft-cd  varied  by  almost  two  orders 
of  magnitude  by  varying  the  amount  of  Cl  doping.  This  feature  allowed  us 
to  use  this  same  type  of  photoconductor  in  assemblies  having  TFEL  output 
layers . 


5.3  Assembly  Techniques 

5,3.1  Thick  Film  EL  Assemblies 


Figure  16  shows  a  cross-sectional  view  of  the  structure 
we  use  for  thick  film  EL  assemblies.  This  structure  is  basically  the 
same  as  we  have  successfully  used  for  some  time.  The  device  is  assembled 
from  the  PC  layer  forward.  The  thick  film  EL  layer  is  simply  sprayed  on. 

5.3*2  TFEL  Assemblies 


This  proved  to  be  a  very  difficult  taek.  We  were  able 
to  dope  arrays  of  CdS,Se  photoconductive  cells  to  give  us  Impedance 
ranges  compatible  with  simple  "sandwich"  electroded  TFEL  lamps.  We  had 
great  difficulty  in  successfully  attaching  our  standard  Interconnect  layer 
to  the  TFEL  lamp  array.  The  interconnect  layer  is  a  glass -epoxy  composite 
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layer  0.003”  thick  which  contains  solid  conductive  posts  on  0.020”  center- 
to-center  spacing  to  connect  the  center  electrode  of  each  PC  cell  to  the 
rear  electrode  of  each  TFEL  cell.  Our  fundamental  problem  seemed  to  be 
that  in  the  initial  turn-on  of  the  TFEL  lamp  array,  there  are  always  a 
number  of  pinhole  defects  that  must  blow  out  and  clear  the  short  circuits 
at  those  points.  Some  of  these  defects  occur  in  the  interface  between 
the  TFEL  lamp  and  the  conductive  post.  It  appears  that  these  defects 
enlarge  to  destroy  the  entire  post-TFEL  interface  and  then  the  entire 
TFEL  cell  breaks  down  and  the  effects  propagate  (probably  because  of  trapped 
gases  and  heat)  to  neighboring  cells.  Eventually,  most  of  the  TFEL  lamp 
array  is  destroyed. 

Our  next  approach  to  the  problem  was  to  omit  the 
interconnect  layer.  We  fabricated  the  PC  cells  on  a  sheet  of  0.003" 
glass  which  had  holes  etched  through  the  glass  in  the  center  of  each 
center  electrode  of  the  PC  cell.  We  bonded  this  PC  substrate  to  the  rear 
of  the  TFEL  lamp  array  with  the  PC  itself  facing  away  from  the  TFEL. 

The  PC  substrate  then  became  the  interconnect  layer.  Then  thin  film 
through-hole  interconnect  techniques  were  used  to  connect  the  PC  cell 
electrodes  to  the  TFEL  lamp  electrodes. 

This  technique  showed  promise.  For  this  technique  to 
work,  an  opaque  dielectric  layer  must  be  included  at  the  rear  of  the  TFEL 
layer  to  prevent  light  feedback  from  reaching  the  PC  layer.  The  best 
thin  film  opaque  dielectric  layer  that  we  made  was  made  of  Ge  co-deposited 
with  Ge02.  l^e  have  achieved  0.1%  transparency  (in  the  green)  with  a 
film  ~  5,OOOA  thick.  The  film  has  a  resistivity  of  greater  than  10^  ohms 
per  square.  The  film  contained  ~  15%  excess  Ge.  We  experienced  problems 
in  achieving  uniformity  of  resistivity  of  the  films.  The  basic  technique 
seems  to  be  that  the  GeOa  evaporation  must  begin  first  with  the  Ge  evapor¬ 
ation  phased  in  slowly  until  the  desired  ratio  of  rates  is  achieved.  The 
usual  results  were  that  high  resistance  films  peeled  off  in  baking  and 
low  resistance  films  adhered  nicely.  The  problem  seems  to  be  in  finding 
the  best  evaporation  rate  schedule  for  getting  both  adhesion  and  high 
resistivity.  We  could  not  complete  this  task  before  the  program  was 
finished.  So  we  fabricated  a  TFEL  light  amplifier  using  a  thick  film 
opaque  dielectric.  This  structure  is  shown  in  Figure  17.  This  structure 
is  cumbersome  but  it  verified  that  the  CdS,Se:CuCl  photoconductor  in 
thin  film  form  can  successfully  control  TFEL  sandwich  cells. 

5.4  Assembly  Test  Results 
5.4.1  TFEL  Assembly 

A  TFEL  assembly  built  as  in  Figure  17  had  response 
characteristics  as  shown  in  Figure  18.  The  points  plotted  are  the 
averaged  output  of  the  assembly.  Because  of  the  sharp  "turn-on  threshold" 
of  the  device,  small  variations  in  TFEL  characteristics  or  PC  characteristics 
cause  the  vertical  portion  of  the  curve  to  vary  over  two-thirds  of  decade 
of  input  values.  We  plotted  the  prediction  vs.  the  average  of  several  cells 
since  the  prediction  is  an  average  of  PC  characteristics  measured  before 
assembly.  The  discrepancy  in  the  dark  output  for  the  assembly  is  most 
probably  due  to  halation  from  nearby  shorted  cells. 
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5.4.2  Thick  Film  EL  Assembl 


A  thick  film  EL  assembly  built  as  in  Figure  16  had  response 
characteristics  as  shown  in  Figure  19.  This  assembly  did  not  fit  the 
computer  simulation.  Investigating  the  discrepancy,  we  discovered  that 
the  opaque  dielectric  layer  in  the  thick  film  EL  layer  had  been  inadvertently 
sprayed  on  too  thick.  This  decreased  the  response  of  the  EL  layer  and 
consequently  diminished  the  gain  we  had  expected.  'rhis  points  up  the 
necessity  for  further  development  of  thin  film  opaque  dielectric  layers. 

The  opaque  layer  should  be  minimized  in  thickness  as  much  as  possible 
so  that  as  much  as  possible  of  the  voltage  applied  to  the  EL  cell  can  be 
across  the  EL  material  and  not  be  wasted  in  the  opaque  layer. 

5.5  Feasibility  of  Stacking  Visible-to-Visible  Image  Conversion 
Panels  onto  Other  Panels 


The  final  assemblies  of  TFEL  and  thick  film  EL  image  conversion 
panels  were  not  finished  in  a  form  compatible  with  actual  stacking  of  one 
onto  another  because  of  electroding  problems. 


Preliminary  simulations  have  been  performed  using  our  best 
CdSe  -  thick  film  EL  combination  from  the  previous  contract  and  data  we 
have  from  previous  visible-to-visible  image  converters.  The  resultant 
combination  would  have  an  output  of  less  than  10 ft-L  for  an  input 
of  ~  8  x  10”6  w/cm2  and  an  output  of  ~  ft-L  for  an  input  of  ~  2  x  10”5 

w/cm2.  The  usable  contrast  range  is  very  narrow  which  makes  the  device 
almost  a  saturating  threshold  type  image  converter.  This  type  of  perfor¬ 
mance  is  due  to  the  fact  that  both  CdSe  and  CdS  PC /EL  image  converters 
have  a  ratio  of  output/input  that  increases  rapidly  to  saturation  as  a 
function  of  input.  When  two  of  these  ratios  are  multiplied  together 
in  a  "stacked"  device,  the  resultant  ratio  yields  and  input-output  curve 
that  resembles  a  step  function.  Figure  20  is  a  plot  of  the  IR  input 
vs .  output  of  a  simulated  device. 


There  are  two  methods  to  overcome  this  lack  of  gray  scale 
range.  The  first  technique  is  to  develop  photoconductors  which  have 
smaller  conductivity  change  per  unit  increase  in  input  radiation  than  the 
present  photoconductors  used.  Past  evidence  indicates  that  some  improvement 
can  be  made  in  this  area. 


The  other  technique  is  to  make  use  of  the  variation  in  input- 
output  properties  of  the  EL/PC  devices  as  frequency  of  voltage  is  changed. 
Figures  21  and  22  illustrate  how  the  input-output  curves  shift  with 
voltage  and  frequency  variation. 

We  are  attempting  to  determine  whether  the  easily  implementable 
techniques  of  amplitude  modulating  and/or  frequency  modulating  the  driving 
voltage  of  the  image  converters  will  effectively  expand  the  gray  scale 
range.  The  basic  idea  is  to  time  share  the  different  input-output  curves 
within  the  limits  of  the  integration  time  of  the  eye  so  that  for  a  fixed 
input  level  a  non-flickering  image  is  perceived. 


Output ,  ft- 


Input,  ft-cd 


PERFORMANCE  OF  THICK  FILM  EL  VI SI BLE-TO- VISIBLE  IMAGE  CONVERSION 

PANEL 


FIGURE  19 


33 


Light  Output,  Ft-L 


Calculated  Overall  IR/Vi3ible  Performance 
of  Two-Stage  Cascaded  Image  Converter 


FIGURE  20 


34 


Light  out,  Ft-L 


10 


1 


Input  Radiation,  w/cm2  (.72^-3. 5^,  2750°K) 


(P  is  the  ratio  of  EL  capacitance  to  shunt 
capacitance  across  the  photoconductor) 

Simulation  of  EL/PC  Panel  Operation  at  100  V  rms 


FIGURE  22 


36 


Toward  this  end,  the  computer  simulations  developed  under  the 
previous  contract  were  "streamlined"  and  expanded  to  assist  in  easier 
reduction  of  greater  quantities  of  data.  The  multi-frequency  driving 
techniques  being  studied  under  a  concurrent  ONR  contract  have  been  simu¬ 
lated  for  an  equal  mixture  of  three  frequencies  and  is  plotted  as  the 
dashed  curve  in  Figure  23.  This  driving  technique  has  been  experimentally 
verified. 


Figure  24  is  the  same  driving  technique  used  with  a  simulation 
of  the  same  CdSe  as  Figure  23  only  using  TFEL  for  the  output  lamps. 

Again  the  three  frequencies  shown  were  mixed  in  equal  proportions  in  the 
composite  simulation. 

It  appears  that  if  the  output  of  the  final  image  conversion 
panel  is  to  be  viewed  and  analyzed  by  the  human  eye,  then  the  slope  of 
the  light-in  vs.  light-out  curve  of  the  final  panel  can  be  effectively  made 
unity.  This  will  avoid  shrinking  the  response  range  of  a  cascaded 
pair  of  image  conversion  panels. 
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6.0  IR-to-Visible  Image  Conversion  Panels 


6.1  Introduction 


Under  this  task  we  were  to  continue  development  of  the  IR 
sensitive  CdSe-thlck  film  EL  panels  we  worked  on  under  the  previous 
contract  (N00014-70-C-0213) .  We  were  to  also  extend  the  work  to  using 
TFEL  as  an  output  media.  Also,  we  were  to  investigate  the  feasibility 
of  placing  a  visible-to-visible  light  amplifier  panel  on  top  of  an  IR- 
to-visible  image  conversion  panel.  In  this  way,  rather  large  conversion 
gains  might  be  realized  for  the  cascaded  pair  of  panels. 

6.2  Results 


We  finished  the  last  program  by  using  a  CdSe  PC  cell  array 
which  was  made  by  depositing  pure  CdSe  onto  an  etched  molybdenum  electrode 
array  (see  Figure  1)  on  a  subgtrate  of  0.003"  thick  glass.  The  CdSe 
was  doped  with  HCl  gas  at  400  C  for  18  minutes. 

Our  first  attempts  at  improving  upon  this  technique  was  to  dope 
the  CdSe  powder  with  CuCl2  solution,  sinter  in  vacuum,  regrind  the  powder, 
and  use  this  powder  as  deposition  material.  The  Cl  would  mostly  escape 
during  sintering  so  the  deposited  films  were  doped  with  HCl  gas  as  before. 
We  were  never  able  to  get  reproducible  or  uniform  films  with  this  process. 
Small  regions  of  some  of  the  films  were  extremely  sensitive  but  uniformity 
was  very  poor.  Vacuum  annealing  of  the  CdSe  films  before  Cl  doping 
didn't  Improve  the  results. 

Next  we  tried  depositing  pure  CdSe,  ~  10, 000 A  thick,  and  then 
depositing  approximately  10 K  of  Cu  on  the  CdSe  films.  The  10 A  of  Cu 
was  deposited  by  an  extrapolation  method.  First  ribbons  of  Molybdenum 
with  a  known  thickness  of  Cu  were  used  as  flash  evaporation  sources  and 
the  thickness  of  the  resulting  deposit  was  measured.  This  measurement 
gave  us  a  source  film  to  substrate  film  ratio  per  unit  area  of  the  source 
ribbon.  From  this  data,  we  prepared  ribbon  sources  to  give  us  a  10 A,  gu 
deposit  on  the  CdSe.  The  CdSe  films  were  then  baked  in  vacuum  at  400  C 
for  1  hour  to  diffuse  in  the  Cu. 

Again,  the  resulting  films  exhibited  very  poor  uniformity. 

Next,  wc  attempted  to  improve  upon  a  procedure  that  was  first 
tried  under  the  previous  contract.  This  procedure  was  to  imbed  the 
CdSe  coated  substrate  in  a  finely  ground  powder  consisting  of  987.  CdSe, 

17.  CuCl2  and  17.  NaCl  as  a  flux.  This  combination  of  imbedding  powder  and 
substrate  was  then  placed  inQa  closed  crucible  and  baked  at  temperatures 
equal  to  or  greater  than  350  C  for  various  periods  of  time.  All  attempts 
produced  films  with  very  non-uniform  response.  We  substituted  Cu2S  for 
the  CuCl2  and  used  various  amounts  of  NaCl  and  many  combinations  of  baking 
times  and  temperatures.  No  improvement  in  uniformity  was  noted.  The 
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literature  had  many  excellent  descriptions  of  techniques  similar  to  this 
being  used  with  good  results.  However  most  experimenters  reported  a 
large  class  of  different  types  of  recrystallization  phenomena  seen  when 
using  this  imbedding  technique.  We  deduced  that  most  previous  experimenters 
had  used  rather  macroscopic  electrode  sizes  and  spacings  when  measuring 
the  response  of  their  resultant  films.  We  have  found  that  we  can  only 
use  electrode  gaps  of  ~  .0015  to  ~  0.004"  in  the  simple  square  "moated" 
electrode  pattern  (see  Figure  1)  and  achieve  proper  ranges  of  cell 
impedance.  The  best  conclusion  we  can  draw  from  all  the  unsuccessful 
doping-by-imbedding  experiments  is  that  this  technique  is  unsatisfactory 
when  using  small  electrode  spacings  and  buried  electrodes. 

We  then  returned  to  trying  by  other  means  to  get  a  CdSe  film 
properly  doped  with  Cu  and  with  Cl  as  a  co-activator.  Knowing  that  using 
a  pre-doped  powder  for  film  deposition  yielded  good  results  for  visible 
light  photoconductors  of  80%  CdS  and  20%  CdSe,  we  again  attempted  to 
duplicate  this  technique  using  only  CdSe  doped  with  Cu  and  Cl.  Initial 
results  from  using  the  pre-doped  powder  were  disappointing.  As  a  point 
of  reference,  thougn,  it  took  more  than  six  man-months  of  dogged  effort 
to  develop  a  reliable  technique  for  depositing  the  pre-doped  CdS,Se  film. 

The  first  few  pre-doped  CdSe  films  exhibited  poor  adhesion  in  patterns 
very  reminiscent  of  the  ones  seen  in  early  attempts  with  the  CdS,  Se  films. 
However,  we  did  not  have  6  man-months  of  effort  left  in  the  program  at 
this  time. 


At  this  time,  we  decided  to  see  just  how  1R  sensitive  our 
green-sensitive  CdS,  Se  films  were.  A  good  green-sensitive  thick 
film  EL  with  CdS,Se:CuCl  assembly  was  tested  on  the  IR  source  and  we 
received  a  pleasant  surprise.  It  was  as  sensitive  as  was  predicted 
for  the  best  cells  of  CdSe:CuCl  that  we  had  occasionally  found  on  some 
substrates.  We  had  never  thought  to  measure  the  IR  response  of  the  PC 
for  this  assembly.  Figure  25  shows  the  response  of  this  assembly  to 
near  IR  and  for  comparison,  the  best  assembly  that  we  fabricated  under 
the  previous  contract. 

Following  this  lucky  turn  of  events,  we  then  tested  a  TFEL  green- 
sensitive  test  assembly  on  the  IR  source.  The  results  are  shown  in  Figure 
26.  Again  the  PC  used  in  this  assembly  had  not  been  characterized  with 
IR  stimulation. 

It  now  appears  that  the  shortest  development  route  to  the 
best  IR-to-visible  image  converter  is  to  modify  our  CdS,  Se:CuCl  deposition 
powder  by  changing  the  ratio  of  CdS  to  CdSe  in  small  increments  until 
the  best  IR  response  is  found. 
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FIGURE  26 


7.0  General  Pfeaults,  Conclusions,  and  Recommendations 

1)  We  expanded  modeling  techniques  to  facilitate  parameterization  of 
materials,  to  investigate  more  complex  forms  of  excitation 
waveforms,  and  to  simulate  the  cascaded  operation  of  various 
potential  PC/EL  candidate  systems. 

2)  TFEL  developed  to  point  where  can  fabricate  in  useful 
brightnesses  (>  50  F-L)  for  useful  lifetimes  (~  1200  hours)  ; 
test  assemblies  have  been  successfully  Integrated  with  photo¬ 
conductors  using  these  films  in  mosaic  form. 

3)  Visible:  P/C  has  been  developed  which,  with  useful  uniformity 
over  large  areas,  is  capable  of  controlling  output  for  input 
values  around  5  x  10"4  foot-candles.  Peak  optical  gains  of  the 
order  of  500  have  been  achieved  using  thick  film  EL  outputs 
(Figure  17),  and  of  the  order  of  3,000  using  thin  film  EL 
outputs  (Figure  18). 

4)  IR  (Figure  25):  Response  curve  moved  down  2  orders  of 
magnitude  with  useful  response  now  available  at  10”6  w/cm2. 

5)  UV:  Unsuccessful  in  obtaining  photoconductors  with  useful 
characteristics,  in  terms  of  impedance  and  voltage  levels 
required  to  control  series  EL  mosaic  cells,  with  measurable 
sensitivity  in  the  UV  spectrum. 

Recommendations : 

Further  studies  of  electroding  and  interconnecting  TFEL  would  be 
required  before  mosaic  arrays  with  the  brightness  and  contrast  required 
for  high-ambient  displays,  such  as  avionic  systems,  can  be  realized  with 
practical  geometries  and  useful  yields. 

State  of  the  art,  in  simulation,  analysis,  design,  and  fabrication, 
has  developed  to  the  point  at  which  specific  applications  in  IR  and  visible 
image  amplifiers,  image  conversion  panels,  and  opto-electronic  logic 
can  be  evaluated  using  simulation  techniques.  In  those  cases  where  a 
simulation  shows  predicted  useful  performance,  fabrication  techniques 
have  been  perfected  to  the  point  where  we  can  predict  the  feasibility 
of  realizing  a  useful  structure. 
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APPENDIX  A 

IR  PROGRAM  AND  PLOT  SUBROUTINE 


A1 .  Introduction 

The  IR  Program  is  written  in  FORTRAN  IV  and  is  used  on  a  time 
sharing  computer  network  with  a  remote  teletypewriter  terminal.  The 
input  files  and  keyboard  inputs  and  the  styles  of  data  output  have 
all  been  programmed  with  ease  of  data  handling  for  the  designer  in 
mind.  The  details  of  the  equations  used  in  the  circuit  analysis  are 
given  in  Section  A2  below  and  a  detailed  description  of  the  program 
itself  is  also  given.  In  addition,  tables  of  all  the  input  parameters 
and  a  'form  of  useful  data  sheets  for  these  input  parameters  are  given. 

A2.  Circuit  Analysis  and  Solution  Procedure 


Figure  Al  shows  the  form  of  circuit  that  this  program  solves. 

It  is  only  a  series  parallel  combination  of  resistors;  however,  the 
resistors  vary  nonlinearly.  We  treat  this  circuit  as  an  almost 
lineal’  circuit.  The  voltage  across  the  EL  cell  is  given  in  rms  form 
by  the  equation: 


v2 


Z1  +  z2 


Equation  Al 


where  V0  is  the  rms  applied  voltage,  is  the  complex  impedance  of 
the  photoconductor,  and  Z2  is  the  complex  impedance  of  the  EL  cell. 
For  simplification: 


let  Sj  -  ou  Cjl  and  S2  **  ui  C2 


then  Z,  -  _ 1 _  -  G1  "  jsl 

G1  4  jSl  Gj2  4  Sx2 

z  .  1  _  G2  -  JS2 

2  G2  +  JS2  G  2  4  S  2 

2  1 

Taking  the  absolute  values,  we  have: 

.1  Zjl  2  -  (Gj2  +  Sj,2)-1 

I  za/  2  “  <°22  ♦ 

I  z  +  Z  i  2  -  (°1  +  +  <S!  *  S2)2 

1  2  COj2  +  Sj2)  (G22  +  S22) 


-  PC  Impedance 

Equation  A2 

-  EL  Impedance 


Equation  A3 
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Circuit  Used  in  Steady  State  EL. 'PC  Image 
Converter  Simulation 
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Putting  these  values  in  Equation  Al,  we  have: 


To  help  in  simplification  in  terms  later  on: 


let  p 


Using  this  simplification,  Equation  Al  becomes: 


Equation  A4 


Equation  A5 


Equation  A6 


Also,  for  the  ratio  of  voltages  of  the  EL  to  the  PC  element,  we  have 


The  photoconductor  conductivity  G^  is  usually  described  by  an 
equation  of  the  form: 

G1  "  u<ud  +  Ls)  Vi*"1  Equation  A9 

Where  U  is  a  dimensional  constant,  is  a  dark  conductivity  term, 

L  is  the  input  radiation,  and  s  and  r  are  constants.  For  complete¬ 
ness,  the  leakage  conductivity  of  the  EL  is  usually  found  to  be  best 
fit  by  an  equation  of  the  form: 

*  H  (wA  exp(b(V2))  Equation  A10 

« 

Where  H  is  a  dimensional  constant,  and  A  and  b  are  constants.  Usually 
b  is  very  small  so  in  the  computations  that  follow,  we  neglect  the 
voltage  dependence  of  Go.  Strictly  speaking,  if  one  were  to  have  a 
large  b  value  one  should  modify  the  derivation  that  follows  to  include 
this  voltage  depcndanco  in  Go.  The  final  thing  to  be  calculated  is 
the  brightness  of  the  EL.  This  brightness  usually  has  a  functional 
dependence  as  given  by: 

n 

Where  D  is  a  dimensional  constant,  F  is  the  driving  frequency,  a  and 
A  are  constants.  Taking  Equations  A9  and  A10  and  substituting  in 
equation  Vj,  we  have: 

G1  “  u(ud  +  LS) (f(G1))r“1  Equation  A12 

Where  f(Gi)  represents  the  right  hand  side  of  Equation  A7 .  In  solving 
this  implicit  equation,  we  now  need  the  services  of  the  computer. 
Equation  All  is  of  the  form: 

Gj  -  f (G^)  Equation  A13 

This  equation  can  be  solved  by  Newton’s  iteration  method. 


) 


Equation  All 


Basically,  this  method  finds  the  zeros  of  a  function  of  the  type: 


y  -  f(x) 

By  repeatedly  doing  an  iteration  of  the  form: 

Equation 

A14 

xn+l  _  xn 

Equation 

A15 

■ —4 


To  get  faster  convergence,  we  used  the  second  order  form  of  this 
iteration  equation  as  given  below: 


y(xn) 

Xn+1  "  xn  “ 


1  +  y(xn)  y '  ’(xn) 
J2  y  ’ (xn))^ 


Equation  A16 


After  the  iteration  is  complete,  the  value  of  Gi  is  used  in  Equation 
A6  to  get  the  rms  voltage  applied  to  the  EL  cell.  The  actual  equations 
used  in  the  second  order  iteration  are  very  long  and  are  contained  in 
the  iteration  calculation  section  of  the  IR  Program. 


A3.  Description  of  IR  Program  and  Plot  Subroutine 


IR  is  written  in  FORTRAN  IV  and  is  used  on  a  time-sharing 
computer  system  through  a  standard  teletype  terminal.  Figure  A2 
contains  a  listing  of  IR.  Figure  A3  contains  a  listing  of  PLOT. 

The  time  sharing  system  it  has  been  used  on  has  excellent  file 
editing  capabilities  if  the  file  contains  line  numbers.  The  first 
field  of  numbers  on  a  line  are  the  line  numbers.  The  second  field 
of  numbers  is  the  FORTRAN  statement  number  field.  The  "*"  are  for 
comment  lines.  The  program  is  liberally  laced  with  comments  as  to 
the  function  of  flags  and  switches  and  the  purpose  of  a  group  of 
statements.  In  the  program,  the  values  of  input  radiation  are 
referred  to  as  "steps".  Unless  problems  with  the  PC  conductivity 
iteration  are  encountered,  the  input  radiation  as  incremented  as 

integers  times  powers  of  ten  i.e.,  1  x  10-6,  2  x  10”6,  3  x  10“6 . . 

1  x  10“5,  .  If,  for  example,  iteration  convergence  is  not 

obtained  for  an  input  step  of  9  x  10“4,  then  the  input  radiation 
value  is  "backed  up"  to  8.1  x  10“4  (8.0  x  10-4  would  have  been 
successfully  used).  If  iteration  convergence  is  still  not  achieved, 
the  input  radiation  is  "backed  up"  to  8.01  x  10“4.  If  convergence 
is  then  achieved,  the  input  steps  are  then  8.02  x  10**4, 

8.09  x  10-4,  8.1  x  10~4,  . .  8.9  x  10"4,  9  x  10“4.  When  the 

program  is  working  in  one  of  these  modes,  the  terms  "back  up"  or 
"small  step"  appear  in  the  comment  lines. 

A  feature  which  might  be  overlooked  upon  first  glance  at  the 
program  or  the  block  diagram  is  that  whatever  form  of  data  output 
is  chosen,  the  program  computes  a  value  of  light  output  for  zero 
radiation  input  for  every  voltage  and  frequency.  This  computation 
is  done  at  the  beginning  of  every  input  light  loop. 

The  linearized  circuit  equations  discussed  in  the  previous 
section  are  utilized  in  the  program  under  the  comment  "iteration 
calculations"  and  at  line  1260  and  line  1300. 

Table  AI  lists  all  the  parameters  for  the  data  files,  the  file 
numbers,  and  gives  a  brief  description  of  the  parameters. 

If  desired,  the  resultant  data  may  be  written  into  a  data  file 
defined  at  program  execution  time. 


Figure  A4  is  a  block  diagram  of  the  IR  program.  This  is  not 
a  flow  chart.  It  is  meant  to  be  used  as  a  guide  to  finding  one's 
way  through  the  program  in  addition  to  providing  a  brief  overview 
of  the  operation  of  the  program. 

Figure  A5  is  a  data  sheet  that  has  been  useful  to  the  computer 
operator  when  using  the  program  and  it  contans  the  format  specifi¬ 
cations  used  in  the  files. 

The  subroutine  PLOT  plots  logarithmically  the  light  output 
versus  the  radiation  input.  The  abscissa  is  not  internally  scaled 
and  accepts  only  the  fifty-five  "normal  step"  values  generated  in 
IR.  The  ordinate  scale  is  controlled  by  the  variable  LP  which 
allows  a  vertical  scale  on  the  teletypewriter  of  six  orders  of 
magnitude.  After  the  plot,  3  lines  of  data  identification  will 
be  printed  if  they  were  included  in  the  data  files.  This  subroutine 
is  used  with  several  other  programs  that  are  described  later. 
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*  IMAGE  INTENSIFIER*  G=U* C UD+T**S>* C V 1 **C R- 1 . 0) > 

*  CHOOSE  TYPE  OF  OUTPUT)  KKK= 1 *  DATA  TABLE 

*  KKK=2*  PLOT  ONLY 

*  KKK=3*  DATA  TABLE  AND  PLOT 
DIMENSION  PC5)*VC5)*FRC5) 

DIMENSION  TT<19*5)*B8CS9*5) 

COMMON  IT*  BB*LP 
DIMENSION  GVOL <  55  ) 

DIMENSION  DVC  3 ) * AV C3)*VEA<3)* ALF  <  3  > 

CALL  DEF IN  E  C  1  >  3HF.L  *  ) 

CALL  DEF I NEC  2* 3HPC*  ) 

CALL  DEF I NEC  3*  6HKN0BS*  > 

CALL  DEF  I  NEC  4#  6M£isI  TE  ) 

KILLD=0 
KKK=  2 

READ  (2*1)  R*  S*  U*  UD 
IFCIE0FC2J.NE. 1 )  READ  <2*552) 

552  FUhMATC •  •> 

KEADC 1*2)  C 

READ  < 1>901)  AALF,  DO*  AA 
IF ( I EOF < 1 >  *NE. 1 )  KEADC 1*551 ) 

551  FORMAT <  *  •> 

READ  (3*4)  CVCJ)*  J= 1  *  5) 

KEADC  3*  4)  (FRC  K ) *  K«l*5) 

READ<  4* 106)  CVEAC I )* 1= 1* 3) 

READC 4* 106)  < DVC I ) *  I  =  1  *  3) 

READC  4* 106)  C  AVC I)* 1  —  1*3) 

READC  4* 1 06)  < ALF C I ) *  I* 1  *  3 ) 

IFC I EOF  C  4) »NE« 1 )  READC4*550> 

550  FORMATC  '  •> 

WRI TEC  9*81 1 ) 

811  FORMATC/ ’CREATE  A  DATA  FILE?  O»N0*  1*YES'/) 
READC  9*812)  KD 

812  FORMATC  ID 
IFCKD.E3.1)  WRITE  <9*S13> 

813  FORMATC/ ’SUPPRESS  UTHER  OUTPUT?  0*N0*  1»YES'/ > 
IFCKD.EQ.l)  READC  9*812)  KILLD 
IFCK1LLD.E0. I >  GO  TO  701 

WRI TEC  9* 107) 

READC 9* 999)  KKK 
705  FORMATC I  10) 

701  tohl TEC  9*  997 ) 

READC 9* 996)  CSHUNT 
CSHUNT*CSHUNT* 1 » 0E~ 1 2 
WRI TEC  9*  5) 

READC 9* 996)  TBEGZN 
LP*AL0G10C  TBEGIN) 

IFCKD.EQ.l)  WRI TE ( 5*  792)  LP 
792  FORMATC I2> 

WRITE  C  9* I  01 )  R*  S*  U* UD 
WRI TEC  9*102)  C 
WRITEC9*  902)  AALF*  DD*  AA 
PC  1 )»C/CSHUNT 
WR 1  TEC  9*3)  PCI) 

WRITE  C  9* 104)  V 
WRITE  C 9# 1 05)  FR 

■  -7 
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1  1  40 
1  1  60 
1180 
1200 
1220 
1240 
1260 
1280 
1300 
1320 
1340 
1360 
1380 
1  400 
1420 
1440 
1  460 
1  480 
lb  00 
1520 
1540 
1560 
1580 
1600 
1620 
1640 
1660 
1680 
1700 
1720 
1740 
1760 
1780 
1800 
1820 
1840 
I860 
1880 
1900 
1920 
1940 
1960 
1980 
2000 
2020 
2040 
2060 
2080 
2100 
2120 
21  40 
2160 
2180 
2200 
2220 
2240 


WRITE(9*517>  DV 
WRITE <9* 518)  AV 
WR1  TEC  9*519)  VEA 
WRITE(9*520)  ALP 

*  NVT-2  PUT  G  VALUES  INTO  GV3LCL)  ARRAY 

*  NVT*1  USE  GVOL(L)  ARRAY  FOR  G  IN  IMPEDANCE  CALC. 
XFCABSCK- 1 • ) »LE. . 1 )  NVT=2 

*  *NCK *  SAVES  LAST  'L*  INDEX  OF  INPUT  LIGHT 

*  LOUP  IN  CASE  OF  SATURATION  BREAKOUT 
t\'CK  =  G 

*  ‘NVTCK  *  USED  AS  SWITCH  TO  STATEMENT  581  IN  CASE 

*  OF  ITERATION  BLOWUP  WHEN  COMPUTING  GVOL(L)  VALUES 

*  USING  OTHER  THAN  FIRST  FKEO*  VALUE 

NV’I CK= 1 

*  SET  '  I VE  *  TO  FIRST  SECTION  OF  B  VS  V  F.L  CURVE 

1  VE=  1 
1  =  1 

*  VOLTAGE  VALUES  DO  LOOP 
DO  88  J  =  1  *  5 

IF  (V(J)-O.O)  99*  99*  10 
10  CONTINUE 

*  PUT  NON-ZERO  VALUES  IN  INPUT*GU1PUT  LIGHT  MATRICES 

DO  516  LM= 1  *  19 
DO  515  LLt)=  1*5 
TT(LM*LLm)«1 .E-7 

515  DB(LM*LLM)s 1 • E-7 

516  CoN’ilNUE 

*  FREQUENCY  VALUES  DO  LOOP 
DO  77  K= 1  *  5 

*  SET  LIGHT  LOOP  TERMINAL  PARAMETER 
I  DK=  1 

IF(FRCK)-0. 0)89*89* 6 
6  CONTINUE 

IF<  KILLD. EQ.  1 )  GO  TO  815 
WRI  TE  (9*  7)  P (  I  ) *  V(U>*  FR(K> 

815  CONTINUE 
301  CONTINUE 

*  PRINT  COLUMN  HEADINGS 

IRC KILLD. EQ» 1 >  GO  TO  922 
I F( I UK*  EO. 1 >  GJ  TO  921 
GO  TO  (921*922*921)*  KKK 

921  WKI TEC  9* 1 2) 

922  CONTINUE 

*  SET  FIRST  LIGHT  VALUE 

T-TBEGIN 

IF(IDK.EQ.l)  T=0. 

*  SET  VALUES  FOR  LOG.  INPUT  LIGHT  VALUE  GENERATOR 
LT«AL0G10(TBEGIN*.99> 

ILT“0 

*  SET  MAIN  PRINT  COUNTER 

N«0 

*  SET  VPC  AND  VI  FUR  ENTERING  LIGHT  LOOP 

VPC*VC J>/4»0 
V1*V(J)-VPC 

*  SET  SUBSCRIPTS  FOR  FIRST  INPUT#JUTPUT  MATRIX  VALUES 
12=0 


Figure  A2 
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2260  KK=  0 

2280  *  INPUT  LIGHT  LOUP 

2300  DO  IS  L=l# IDK 

2320  *  INPUT  LIGHT  LOG  STEP  GENERATOR 

2340  ILT=ILT+1 

2360  G U  TO  <53#53#53#53#53#53#53#53#53#52)#ILT 

2380  52  LT=LT+1 

2400  ILT=  1 

2420  S3  CONTINUE 

2440  * 

2460  *  SET  ITERATION  BLOWUP  COUNTER 

2480  MM=  1 

2500  *  ADVANCE  NORMAL  LARGE  LIGHT  STEP  PRINT  COUNTER 

2520  N=N+1 

2540  *  SET  SMALL  STEP  PRINT  COUNTER 

2560  MMM= 1 

2580  *  REENTER  FROM  SUCCESSFUL  SMALL  STEP  CALC. 

2600  51  CONTINUE 

2620  *  SET  ITERATION  COUNTER 

2640  M=0 

2660  *  EN 1EK  ITERATION  LOOP  ON  SMALL  LIGHT  STEP 

2680  87  CONTINUE 

2700  IF  <L“ 1 >  23#  23#  22 

2720  22  VPC=V<J>-VE 

2740  V 1  *  VE 


2760  23  CONTINUE 

2780  *  CHECK  VOLTAGE  DEPENDENCE  OF  CONDUCT  I VI  TV 

2800  IFC ABS( K- 1 • ) » GT» • l >  GJ  TO  863 

2820  GO  TO  <  867#  863)#  NVT 

2840  867  IF(L.GT.NCK)  GO  TO  863 

2860  G=  GVUL ( L  > 

2880  GO  TO  864 

2900  863  G=U+(  UD+T**S>*(  VPC**CR-  I  *  >  > 

2920  864  CONTINUE 


2940 

2960 

2980 

3000 

3020 

3040 

3060 

3080 

3100 

3120 

3140 

3160 

3180 

3200 

3220 

3240 

3260 

3280 

3300 

3320 

3340 

3360 


4 

♦ 


*  NORMAL  ITERATION  REENTRY 
8  CONTINUE 

*  COMPUTE  IMPEDANCE  OF  PC 
Y1=P< I)*G/<6.28+FK<K>*C> 

'  ADVANCE  HE  RATION  COUNTER 
M*M+  1 

*  COMPUTE  1/K  OF  EL 

Y2USSDD*(FR(K )**AALF>*EXP(AA*(V<J)-V1 > ) 

*  COMPUTE  POWER  FAC  1  OR 
PF»  1  ./S3RTC  C  (  (  6.  2S*FRCK)*C>/Y2G>**2>M  »  > 

ITERATION  CALCULATIONS 
Y2«P< I >+PK/SORT< 1 • -PF*PF  > 

Y-Y1+Y2 

Z*( 1 .+P< I ) )*( 1 . +P( I >  > 

V1«V< J)4S0KT<<Y2*Y2*P<I>*P(I >)/cy*Y*Z>> 

G1*U*(UD  +  T**S)*CV1 <  R“ | * } ) 

DIFGMK-1 . >*P< I >*G1 / ( 6* 28*KR<K>  *C*<  Y*Y  +  Z) ) 

DDIFG»< ( < DIFG*P< I )/6. 28*FK(K )*C*Y) ) ♦( 2.*C-DIFG)*P( I )*Y* 
CD1FG)  >/C6.  28*FR<K)*C*<Y*Y«-Z>  >  >/<  1  .♦<DIFG/G1  )  ) 
UD«<Gl-G>/<-DIFG-l.) 

G2»G-CGD>*< 1 •♦<CGD)*DDIFG>/<2.*(-DIFG-l . >>> 
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3380 

* 

CONVERGENCE  TEST 

3400 

ESC= 

•ARS(  (G2-G)/G2) 

3420 

GO  TO  (51 1,512,513), 

MM 

3440 

51  1 

1 F ( ESC- 0 «  0 1 )  20, 

20, 

13 

3460 

512 

IFC  ESC-O*  02)  20, 

20, 

13 

3480 

513 

1  F(  ESC-O*  05)  20, 

20, 

13 

3500 

* 

VALUE  FOR  NEXT  ITERATION 

3520 

13 

G=G2 

3540  If  (M  -  500  >  66,  66,  44 

3560  66  GO  TO  8 

3580  *  END  OK  ITERATION,  NO  CONVERGENCE 

3600  *  TRY  ADVANCING  LIGHT  INPUT  BY  1/10  PREVIOUS  STEP 

3620  A  A  GO  TO  C  A3*  131, 132), MM 
36 40  +•  FIRST  BACKUP,  i'iM=  1 
3660  43  T-'i-<(10.0**LT>*0.9> 

3680  M=0 

3700  MM »2 

3720  MMM=2 

3740  NVT=  ft 

3760  GO  TO  87 

3780  *  SECOND  BACKUP,  MM =2 

3800  131  T=  I’-  <  C  10»0**LT>*0»09) 

3820  M=0 

3840  MM =3 

3860  NNN=MMM 

3880  MM  Ms-  2 

3900  GO  TO  37 

3920  *  ITERATION  HAS  CONVERGED 

3940  20  CONTINUE 

3960  *  COMPUTE  OUTPUT  VARIABLES 

3980  Y3=  P (  I  )  * G2/  (  6  .  28 *F  R (  K  >  *C  > 

4000  VE“  V(  J  >  *.SORT(  C  Y3*Y3+1 . 0  > /  (  (  Y3  +  Y2 )  *  C  Y3  +  Y2)  +?.  >  ) 

4010  DO  181  IT*  1,3 

4020  IF(VE-VEACIl)  >  182,181,181 

4030  182  IVEs II 

4035  GO  TO  183 

4040  181  CONTINUE 


4050 

4060 

4080 

4100 

4120 

4130 

4140 

4180 

4200 

4220 

4240 

4260 

4280 

4300 

4320 

4340 

4300 

4380 

4400 


183  CONTINUE 
400  DVLOG=ALOG(  L)VC  I  VE  >  > 

FLJU*ALf  < I VE  >  *( ALOG(FR(K) >  > 

ELOt---SJPTCAV<  IV£)/VE) 

PLOGs  DVL  UG+FL  JG- EL  0  G 
IF C CLOG. LT. -23.026)  BLOG* -23 .026 
B=EXPC  BLOG) 

♦  SMALL  STEP  PRINT  CHECK 

GJ  TO  (393,396,395,396,395,396,395,396,395,394),  MMM 

396  V>KI TE( 9,  11)  T,  B,  VE,  G2 
395  CONTINUE 

GO  TO  (393,397,  398),  MM  - — - 

398  ’1  ■  T *  ( ( 1 0 •  0**L T )  *0 •  0 1 )  Figure  A2 

GO  TO  399 

397  T*T+(  ( 10»0**LT>*0. 1 )  1'a^e  4  of  IR  Program 

399  MMM* MMM ♦ 1  — 

M»0 

GO  TO  87 

394  GO  TO  (571, 572, 573), MM 
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4420 

4440 

4460 

4480 

4S00 

4520 

4540 

4560 

4580 

4600 

4620 

4640 

4660 

4680 

4700 

4720 

4740 

4760 

4780 

4800 

4820 

4840 

4860 

4880 

4900 

4920 

4940 

4960 

4980 

5000 

5020 

5040 

5060 

5080 

5100 

5120 

5140 

5160 

5180 

5200 

5220 

5240 

5260 

5280 

5300 

5320 

5340 

5360 

S380 

5400 

5420 

5440 

5460 

5480 

5500 

5520 


573  T*T+<< 10.**LT>*.01 > 

MMM=NNN 
MM*  2 

GO  T3  571 

572  T=T+<<10.**LT)*.l) 

MMM=1 
MM=  1 

GO  TO  <571/581 )/NVTCK 
581  NVT=J 
571  GO  TO  51 

*  END  OF  SMALL  STEP  CALC. 

*  MAIN  STEP  PRINT  CHECK 
393  CONTINUE 

IK < IDK.EQ* 1  )G0  TO  951 
IF < ABS<  R- 1 »  > • GT* • 1 >  GO  TO  862 
GO  TO  <866/ 869)/NVT 
866  IF<L.LE.NCK>  GO  TO  862 
869  GVOL  <L )  =  G2 

862  GO  TO  <96/ 96/ 15/ 15/ 96/ 1 5/ 1 5/ 1 5/98)/N 
96  CONTINUE 

GO  TO  <693/930/930)/  KKK 
930  CONTINUE 
11*11+1 
JJ=K 
KK=KK+ 1 
LL=K 

TT< II/ JJ)*T 
BB<KK/LL >=B 

IF  <KD*E0»  1  )  WRITE  <5/810)  B/  T#  V<  J)/FR<K) 

810  FORMAT  <4E  12.  4) 

IF< I I- 1 9 )  693/692/692 
692  11*0 
KK*0 

693  CONTINUE 

GO  TO  <951/950/951)/  KKK 
951  CONTINUE 

IF<KILLD.EO.  1 )  GO  TO  814 
WRI TE< 9/11)  T/  B/  VE/  G2 
814  CONTINUE 

IF  < IDK.EQ. 1)  GO  TO  300 
950  CONTINUE 

*  BREAKOUT  OF  LIGHT  LOOP  IF  OUTPUT  IS  SATURATED 

IF  <L-  1  >  232/  232/  231 

231  IF  <<VE/V<J>>-0.995)  232/  111/  111 

232  CONTINUE 

98  IF  < N “ 9 )  612/  61 1/  61 1 

61 1  N*0 

612  CONTINUE 

*  END  OF  LIGHT  STEP  LOOP 
15  T*T  +  < 10.044LT) 

IDK*  1 
NCK-55 
GO  TO  111 

*  CHANGE  LIGHT  LOOP  TO  FULL  RANGE 
300  IDK* 55 

GO  TO  301 
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5540 
5560 
5580 
5600 
5620 
5640 
5660 
5680 
5700 
5720 
5740 
5760 
5780 
5800 
5820 
5840 
5860 
5880 
5900 
5920 
5940 
5960 
5980 
6000 
6020 
6040 
6060 
6080 
6100 
6120 
6140 
6160 
6180 
6200 
6220 
6240 
6260 
6280 
6300 
6320 
6340 
6360 
6380 
6400 
6420 
6440 
6460 
6480 
6500 
6520 
6540 
6560 
6580 
6600 
6620 
66  40 


*  ESCAPE  FROM  LIGHT  LOOP  <  SATURATED  OUTPUT  ) 

111  CONTINUE 

NVT=  1 
NVTCK® 2 

IFCNCK.EG. 55)  GO  TO  63 
NCK=L 

63  CONTINUE 
GO  TO  77 

*  NO  ITERATION  CONVERGENCE  FOR  THIS  FREQ. *  VOLTAGE 
13?  LRI TEC  9*33)  ESC 

*  END  OF  FREQ  DO  LOOP 
77  CONTINUE 

89  CONTINUE 

GO  TO  <940# 941*941 )>  KKK 
941  CONTINUE 

*  TALL  PLOTTING  SUBR. 

IF CKILLD.EQ. 1 >  GO  TO  940 
CALL  LUGPL 
WRITE? 9#  5  52) 

LRI TEC  9* 55.  ) 

WHITE? 9*550) 

940  CONTINUE 
NVT-2 
NCK=0 
NVT  CK= 1 

*  END  OF  VOLT  DO  LOOP 
88  CONTINUE 

99  CONTINUE 

* 

*  FORMAT  STATEMENTS 

517  FORMAT < ' DV  =  •> 3< 1 PE8. 2* 2X> > 

12  FORMAT  ?  'WATT  S/CM2  '*  ?X*  ’LIGHT  0UT’*4X,  ’V  EL  ’  #  6X*  ’  PC0ND  ’  ) 
11  FORMAT? 1P4E10. 2)  . 

1  FORMAT? 4E 15. 2) 

101  FORMAT?// 4K  K»* 1  PE  1 0*  2*  4H  S=*E10.2*4H  U=#1PE10.2* 

+  4H  I.IDa #  1  PE  10*  2) 

2  FORMAT ? El  5* 6) 

102  FORMAT?’  C<='#  1PE10.2) 

901  FORMAT ?3E1 5.6) 

902  FORMAT ?5HAALK«#  1PE9.2/4H  DD=*  1  PE  1  0. 2*  4H  AAs#E10.2> 

3  FORMAT? ’CAP  RAT  1 0= • > 1 PE6. 2) 

518  FORMAT?  'AV  «  '*  3?  1PE8.2*2X)  ) 

4  FORMAT? 5F 10. I ) 

104  FORMAT  ?  1 0)1  V  APPL=  *5F10.!> 

5  FORMAT? ’BEGINNING  INPUT  KAD.'/> 

105  FORMAT?  10H  Fi'.EC  =  *5F10*1> 

7  FORMAT? 10HCAP  KA T IU  =  # F6. 1 » I  OH  V  APPL«/F6. 1# 

♦  OH  F  RED®  *  F7 • 1 ) 

33  FORMAT? 3HESC# 1PE20. 8) 

106  format? 3ei 6. 2) 

107  FORMAT? ’TYPE  ”2"  IF  PLOT  IS  DESIRED*/ 

+  ’TYPE  "l"  IF  ONLY  DATA  IS  DESIRED’/ 

♦  *1 YPE  ”3"  IF  ROTH  DATA  AND  PLOT  ARE  DESIRED’//) 

997  FORMAT? ’TYPE  PC  SHUNT  CAPACITANCE  IN  PF’//> 

996  FOhMATCEl 6. 2) 

999  FORMAT? 1 1 > 


FIGURE  A2 
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XR  PROGRAM 

6660  519  F0RMATC ’VEA  *  % 3< 1PE8. 2, 2X> ) 
6680  520  F0RMAT< 'ALF=  ** 3< 1PE8. 2» 2X>/> 
6780  STOP 

6800  END 
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PLOT  SUBROUTINE 

20  SUBROUTINE  LUGPL 

40  .  INTEGER  X/  PLOT  /  BLANK/  DOT 

'60  DIMENSION  LINE(6l>/  XX<19/5>/  YY  <1  9/  5  )  t  Y  (  5  )  *  PLOT  C  5  > 

80  COMMON  XX/  YY/LP 

100  *  CHOOSE  GLUTTING  SYMBOLS 
120  DATA  BLANK/ 1 H  /t DO  1/ 1 H*// X/ 1  MX/ 

1/0  DATA  PLO  T  /  1 H  1  .*  1  H2/  1  H3/  1  HA*  1  H5/ 

160  *  LABEL  THE  PLOT 
160  kiUTE  C9/P01) 

200  201  FmvMA  1  C///20H  LIGHT  OUT '/  FT  L  > 

220  207  V.fil  TE  C  9*  600) 

2/0  600  FORMAT!  •  1  0-4S  5X/  *  10-3'»6X>  *  1  0-2  '  t  6X/  '10-1  *»  7  At 
260  +  '1  SOX#  *10+1  '/6X/  *10  +  2') 

280  *  PRINT  LEFT  BORDER 
300  99  CONTINUE 

320  DU  101  N- 1,61 

3/0  L  I N  E  (  N  )  =  D  U  ] 

360  101  CONTINUE 
380  602  DO  60S  1-1.7 
ZOO  Lifts 1+10*< 1-1 > 

420  LtNE!LIN)=X 

440  60S  CONTINUE 
460  603  l.  RITE  <9*  102)  LINE/  LP 

480  102  FORMAT  < 6 1  A 1  *  '  1  * OE ' / I  2 )  ! 

SOO  *  BLANK  THE  LINE 
S20  DO  103  1=1/61 

840  LINE! I ) -BLANK 

860  103  CONTINUE 

860  *  kb.  ONLY  PLOT  1  /  2/ 5/ 1  0/ ETC  ON  X-AXIS*  2=  2/ 5=4  /  10=6 
600  *  COMPUTE  YY  POINTS 
620  DO  105  K">  1  /  19 

640  DU  104  L= 1  /  8 

660  IF!  YYCK/L  > . LT» 1 •  E-  7>  YY < K/ L ) = 1 . E- 7 

680  Y ( L )  =  AL 0  G  1  0  (  Y  Y  (  K  /  L  )  ) 

700  222  YY  <  K/  L  )  =  Y  (!..  )  F  1  0  •  +  4 1  • 

720  104  CONTINUE 
740  108  CONTINUE 
760  *  X  AXIS  STEP  LOUP 
730  DU  107  1-1/19 

800  *  MAKE  SURE  LINECN)  IS  BLANK 
820  DO  120  NN  = 1/61 

840  L INE(NN)= BLANK 

860  120  CONTINUE 

830  +  LOAD  IN  PLOTTING  SYMBOLS  IN  CORRECT  POSITION 
900  I CMECK=0 

920  DU  106  J= l  /  8 

940  N=YY(1/J) 

960  IF  <N-0)  106/  106/  103 

980  103  IF  (  N  -  6  1  )  109.  109/  106 

1000  109  LINE(N)=PLJT < J) 

1020  I  CHECK® I  CHECK + 1 

1040  106  CONTINUE 
1060  L  INF!  C  1  )  =  DOT 

1080  *  SELECT  HOKIE.  POSITION 
1100  GO  TO  C  1  0  7  /  1  1  1  /  1  1  2/  1  1  4/  l  I  I  /  1  1  2/  l  1  4/  1  1  1  #  1  l  2*  1  I  4/  I  1  1  /  I  1  2/ 

1120  +  114/111/112/114*111/112/114)/  I 


FIGURE  A3 
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PAGE 

PLOT 

1140 
1160 
1180 
1200 
1220 
1240 
1260 
1280 
1300 
1320 
1340 
1360 
1380 
1400 
1420 
1440 
1460 
1480 
1  500 
1  520 
1540 
1560 
1580 
1600 
1620 
1640 
1660 
1680 
1690 
1700 
1720 
1740 
1760 
1780 
1800 
1820 
1840 
1860 
1880 
1 B90 
1900 
1920 
1  V  40 
1960 
1970 
1980 
2000 


2 

SUBKUUT  1NE 

*  PLOT  2 

111  IF  C I  CHECK-0 )  400*  400*  411 
400  WRITE ( 9#  402)  LINE  C 1 ) 

402  F0RMAT< 1H./1A1 ) 

GO  TO  107 

411  IC-61 

DO  70  I  V=  1  *  6 1 

IF(L1NE( IC) «NE« BLANK)  GO  TO  80 
IC=IC- 1 
70  CONTINUE 
80  CONTINUE 

503  WKITEC9*  403)  (LINE (NNN ) *  NNN=1*IC) 
60  TO  107 

403  F 0 KMAT ( 1 H • / 6 1  A 1 ) 

*  PLOT  5 

112  IF( I CMECK-0 )  400*  400*  412 

412  IC-61 

DO  82  1 V*  1  *  61 

IF(LINE( 1C) .NE. BLANK)  GO  TO  84 
IC=IC-  1 
82  CONTINUE 
84  CONTINUE 

WRITE<  9* 404)  (LINE (NNN >* NNN* 1 » IC> 

404  FOKMAT( 1H./61A1 ) 

GO  TO  107 

*  PLOT  1*  10*  ETC. 

1 1 4  L I N  E ( 1 )  =  X 

LINE( 61 >  =  X 
LP*LP+1 

IF  (1-19)  443*  445#  445 
443  CONTINUE 

IF ( I  CHECK-0 )  400#  400*  444 

*  RIGHT  BORDER  DECADE  MARKERS 


445 

DO  714  I UUM* 1 #  7 

LIN* 1 + 1 0*( I DUM- 1 > 

LINE(LIN)=X 

714 

CONTINUE 

713 

WRITE (9* 301)  LINE* 

LP 

301 

FORMA T( 1H./61 Al#  •  1 

•  OE 

GO  TO  107 

444 

CONTINUE 

51  1 

WRITE  ( 9#  4C5)  LINE 

#  LP 

405 

F  u  KM  A  T  (  1  H .  /  6  1 A 1  *  *  1 

•  OE 

107 

CONTINUE 

LP=LP-6 

RETURN 

END 

FIGURE  A3 
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TABLE  A I 


File 

Name 

EL 


PC 


KNOBS 


Input  Parameters  for  IR  Program 


Input 

Device 

1 


2 


3 


Variable 

C 

AALF 


DD 


AA 


(Optional ) 
R 

S 

U 

UD 


(Optional) 

V(I) 

FR(I) 


_ Comments _ 

Capacitance  of  EL/cell 
Frequency  dependence  of  EL  leakage 
conductivity.  GEL  oc  (Freq.)  AALF 
Dimensional  constant  (H)  for  EL  leakage 
conductivity  (EL  conductivity  is  in 
mhos) 

Coefficient  (b)  of  voltage  dependence  of 
EL  leakage  conductivity.  G£L  oc 
Exp(AA* (Voltage) ) 

30  character  line  for  data  identification 
Exponent  of  voltage  dependence  of  PC 
conductivity.  Gp^oC  (Volt)11”1 
Exponent  of  input  energy  dependence  of 
PC  conductivity.  Gpc  oc  (L)s 
Dimensional  constant  for  PC  conductivity 
(PC  conductivity  in  mhos) 

Dark  current  constant  (U^)  for  PC 
conductivity 

30  character  line  for  data  identification 

Array  of  up  to  5  r.m.s.  applied  voltages 
[in  Volts] 

Array  of  up  to  5  applied  frequencies 
[in  Iiz] 
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TABLE  A I 

Input  Parameters  for  IR  Program  (cont.) 


File  Input 

Name  Device 

BRITE  4 


Keyboard 

Inputs 


Variable 

VEA(I) 


DV(I) 


AV(I) 


ALF(I) 


(Optional ) 
CSHUNT 

TBEGIN 


_ _ Comments _ 

If  brightness  of  EL  vs.  voltage  fits  the 
equation  Bocexp(-  only  over 

limited  ranges  one  can  fit  this 
equation  to  the  experimental  curve 
over  3  ranges  of  voltage.  VEA(I) 
are  the  upper  voltage  limits  of  each 
approximation . 

Array  of  up  to  3  dimensional  constants 
to  fit  an  EL  brightness  vs.  voltage 
curve  (Brightness  in  Ft-L) . 

Array  of  up  to  3  exponential  constants 

to  fit  the  equation,  Boc  Exp(-n  /AVQ)^ 

V  Volt 

to  an  EL  brightness  vs.  voltage  plot 
over  the  full  range  of  possible  applied 
voltages . 

Array  of  up  to  3  exponential  constants  of 
EL  brightness.  B  oc  (Freq.  )^L^^^derived 
from  fitting  of  an  EL  brightness  vs. 
voltage  curve. 

30  character  line  for  data  identification. 

Shunting  capacitance  around  the  PC  due 
to  structure  of  device  [in  pf].  (E16.2) 

Smallest  input  radiation  value  for  the 
simulation.  (£16.2) 
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Inputs : 


EL  Parameters  PC  Parameters  Structural  Cap.  Power  Supp 


BLOCK  DIAGRAM  OF 


IR 


PROGRA 


\t 


FIGURE  A4 


© 


|  "Computo' Concluc't'ivity  of  EL  1 

. — . . H  ■■  . .  -  — 

|  Compute  Power  Factor 


I  Second  Order  Newton  Iteration  Calculation  j 


Convergence  Check? 
(500  Iterations) 


no ' 


yes 


Set  New  PC 
Conductivity 


no 

I  • 

second  500 
attempts 


‘no 

I 

third  500 
attempts 


int  Convergenc 
Difference 


Compute  Output 
Variables 


Decrease  Input  By 
9/10  Previous 
Attempted  Increment 


Try  Another 
Frequency _ 


Again  Decrease  Input 
by  9/10  Previous  Attempted 
Small  Increment  _  _ 


Relax  Criteria 
for  Convergenc o 


3L 


Further  Relax 
Conver gence  Criteria 


1 - > 

* 

_ _  l 

J - - 

Print  Check 
for  Size  of 
Input  Steps 

Set  Flags  and 
Switches  for  Print¬ 
outs  of  Small  Input 
Steps 

Set  Flags  and  Switches 
for  Printouts  of  Smallest 
Input  Steps 

X 

— 

Am  at  Normal 
Large  Step 
Value 


© 


Am  at  a  Small 
Input  Step  Value 


31 


BLOCK  DIAGRAM  OF  IR  PROGRAM 
FIGURE  A4 
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Less  than  Max.  Max. 


BLOCK  DIAGRAM  OF  IR  PROGRAM 
FIGURE  A4 
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INPUT  DATA  FORMAT  FOR  IR 


DATE: 


/ 


EQUATIONS: 

B 


DV(I)  *  FALF(I)  exp  ~ (J) 


Y2G(Leakage)  -  DD*FR(k)AALF  ov„(AA*V( J) -VI) ) 

exp 

G1 (g)  =  U*(UD  +  TS)*(V1R"1) 


DATA 

BLOCKS : 

NAME 

DEVICE 

DATA 

EL 

1:  C 

(E15.6) 

AALF 

DD 

AA 

(E15.6) 

(Optional : 

30 

character 

line 

for 

labeling) 

PC 

2: 

R 

S _ 

U 

UD  (E15.2) 

(Optional : 

30 

character 

line 

for 

labeling) 

KNOBS 

3: 

V(I) 

(F10.1) 

FR(I) 

(F10.1) 

BRITE 

4: 

VEA(I)  _ 

(E16.2) 

DV(I)  _ 

(E16.2) 

AV(I)  _ 

(E16.2) 

ALF(I)  _ 

(E16.2) 

(Optional: 

30 

character 

line 

for 

labeling) 

Test  Data: 


FIGURE  A5 
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APPENDIX  B 


IMINT  PROGRAM  AND  RUNGE  SUBROUTINE 


Bl.  Introduction 


The  IMINT  Program  treats  essentially  the  same  circuit  as  the 
IR  Program  with  one  important  difference.  The  difference  is  that 
the  conductance  of  the  photoconductor  is  allowed  to  vary  exponentially 
with  time  in  response  to  step  functions  light  input.  A  schematic 
diagram  of  the  circuit  treated  is  shown  in  Figure  Bl . 

B2 .  Circuit  Analysis  and  Solution  Procedure 


To  solve  this  circuit  we  do  not  use  any  linear  approximations. 
Instead,  we  write  the  equations  for  the  total  current  through  the 
combined  circuit  in  a  form: 


I 


'EL 


d(VEL) 


dt 


+  VEL 


1 

REL 


Equation  Bl 


I  -  C 


PC 


d(VDP) 


P£i  +  V 


dt 


PC 


R 


PC 


Equation  B2 


The  conductivity  of  the  EL  cell  is  written  again  as  it  was  for  the 
IR  Program  and  is  expressed  below: 


H  FA  eXp(K(VEL)) 


Equation  B3 


whore  H  is  a  dimensional  constant,  F  is  the  frequency,  and  A  and  K 
are  constants.  The  photoconductor  conductivity  is  again  expressed 
in  the  same  form  as  before,  but  modified  with  an  exponential  time 
term. 


In  the  equation  below,  the  upper  exponential  term  is  for  the 
duration  of  the  light  pulse.  The  lower  time  term  is  for  the  decay 
time  of  the  photoconductor  after  the  input  radiation  pulse  has  ceased. 


1 

HpC 


U(Ud  +  Ls) 


r-1  i"  1  -  oxg(-X^(t-t^))f 


if 

if 


t^ctct 

t2<t 


Equation  B4 
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•  V. 


Light 

Out 


/f:  TL exp(bVEL) 


Optical  Barrier  tJ)  Prevent  Feedback 


vRadiation  Input 
(Square  Impulse,  peak  Value 


Before  Radiation  Pulse:  (t<tj) 


Hpc  '  <«><t'd>(Vpc> 

During  Radiation  Pulse:  (t^<t<tg) 

R  -U<0d  +  LB)(Vpcf"l[l-cxp(-/.1(t-t  ))3 

PC  A 

After  Radiation  Pulse:  (t2*:t) 

r  "  u(ud  +  LS)(vpc)r  1Coxp(-j.2(t-t2)) 3 
PC 

Circuit  Used  in  Simulation  of  EL'PC  Image  Converter 
RosponKC  to  Transient  St  inula ti on 


FIGURE  Bl 


•  .0/1 


In  this  equation,  U  is  dimensional  constant,  Ud  is  a  dark  conductivity 
term,  L  is  peak  input  radiation,  s  and  r  are  constants,  Xi  is  a  rise 
time  constant,  X2  *s  a  decay  time  constant,  t^  is  the  beginning  of  the 
light  pulse,  and  t2  is  the  end  of  the  light  pulse.  The  applied 
voltage  is  assumed  to  be  of  the  form: 

Vapp  -  Vm  sin(cut-60)  Equation  B5 

Where  6q  is  chosen  so  that  the  beginning  voltage  on  the  EL  for  the 
simulation  is  zero.  This  gives  us  for  the  voltage  on  the  photocon¬ 
ductor  the  following  term: 

• 

Vpc  -  tv.  sin(tut-60)  3-Vjjl  Equation  B6 


Putting  Equation  B6  into  the  combined  Equation  B1  and  B2,  we  have 
Equation  B7. 

<CEL  +  Cpc)  d  VEL  «  -VEL[HFA  exp(K(VEL))] 
at 

+CpcCVina)  cos(tut-60)  ] 


+  [U(Ud  +  Ls)  ]  [Vm  sin(cut-60)  -  V, 


EL 


l  expC-XgCt-tg) )  J 


Equation  B7 


This  equation  is  of  tho  functional  form: 


Equation  BS 


This  equation  is  a  first  order  equation  which  is  easily  integrable 
by  several  numerical  techniques.  We  have  used  one  of  the  simplest 
techniques  which  is  the  Rungo-Kutta  mothod.  The  present  IMINT 
Program  is  not  ideally  suited  for  tho  simulation  of  thin  film  EL 
light  output  because  it  has  no  provision  yet  for  tho  inclusion  of 
multiple  curve  fits  to  the  TFEL  brightness  versus  voltage  curve. 

B3 .  Description  of  IMINT  Program  and  Itunge  Subroutine 

IMINT  is  written  in  FORTRAN  IV  and  is  used  on  a  time  sharing 
computer  systom  through  a  standard  telotype  terminal.  Figure  B2 
contains  a  listing  of  IMINT.  The  time  sharing  systom  it  has  been 
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used  on  has  excellent  file  editing  capabilities  if  the  file  contains 
line  numbers.  The  first  field  of  numbers  on  a  line  are  the  line 
numbers.  The  second  field  of  numbers  is  the  FORTRAN  statement  number 
field.  The  are  for  comment  lines.  The  program  is  liberally 
laced  with  comments  as  to  the  function  of  flags  and  switches  and 
the  purpose  of  a  group  of  statements.  A  complete  listing  of  the  input 
parameters  for  the  program  along  with  the  device  number  called  for 
in  the  program  are  shown  in  Table  Bl. Table  BII  contains  the  input 
parameters  which  are  entered  at  the  keyboard  at  program  execution 
time  along  with  comments  as  to  the  nature  of  the  parameters.  The 
IMINT  Program  integrates  Equation  B4  in  the  section  above  by  means, 
of  a  subroutine  called  RUNGE.  The  program  is  basically  a  simple 
integration  program  except  for  two  parts.  The  first  of  these  parts 
is  the  determination  of  the  phase  difference  between  the  applied 
voltage  and  the  voltage  on  the  EL.  This  phase  difference  (60)  is 
determined  such  that  the  voltage  on  the  EL  at  the  beginning  of  the 
simulated  voltage  cycle  is  passing  through  zero  in  the  positive 
diref tion.  The  program  finds  this  phase  angle  from  an  initial  guess 
by  the  designer  by  moans  of  an  iteration  technique  and  computing  the 
rate  of  change  of  the  convergence  difference  between  two  cycles  run 
with  different  phase  angles  and  computing  a  new  phase  angle  based  on 
this  rate  of  change  of  convergence  with  respect  to  phase  angle.  Once 
the  program  has  determined  a  60  so  that  the  voltage  on  the  EL  returns 
to  zero  in  one  period  of  the  applied  frequency,  the  program  then 
integrates  from  the  end  of  this  one  cycle  (with  no  input  radiation 
during  this  cycle)  forward  in  time  for  as  many  cycles  as  the  designer 
choses . 

The  second  part  of  the  program  differs  from  a  simple  integration 
routine  because  of  the  time  varying  nature  of  the  radiation  input. 

This  time  varying  nature  is  handled  by  introducing  the  exponential 
rise  of  decay  functions  as  shown  in  Equation  B4  of  the  section  above 
at  the  appropriate  time  in  the  integration. 

Figure  B3  is  a  block  diagram  of  the  IMINT  Program.  This  is  not 
a  flow  chart.  It  is  meant  to  be  used  as  a  guide  for  finding  one's 
way  through  the  program  in  addition  to  providing  a  brief  overview  of 
the  operation  of  the  program.  Table  BIX!  and  Table  BIV  are  data 
sheets  that  have  been  useful  to  the  computer  operator  when  using  the 
program  and  contain  the  format  specifications  used  in  the  files. 

The  subroutine  RUNGE  is  a  simple  form  of  the  Rungo-Kutta  technique 
for  integrating  first  order  differential  equations.  The  flags, 
switches,  and  variables  in  the  subroutine  arc  described  in  the  comment 
statements  at  the  beginning  of  the  subroutine  listing.  A  listing 
of  RUNGE  is  contained  in  Figure  B4, 


inipH^rpaGfipfiv m  i 1 ,i  mm f  i ■  m  mmm  upu*  «  **  *»v 


As  yet,  we  do  not  have  a  useful  plotting  subroutine  to  plot 
the  output  of  the  IMINT  Program.  At  present,  the  output  is  the 
time,  voltage  on  the  EL,  EL  current,  and  derivitive  of  the  voltage 
on  the  EL.  The  data  must  then  be  plotted  by  hand.  From  the  plot, 
one  can  find  peak  AC  voltages  applied  to  the  EL  as  a  function  of 
time  and  use  this  in  a  simple  calculation  with  EL  brightness  formulas 
that  are  appropriate  for  thick  film  or  thin  film  EL  to  get  the 
brightness  output . 
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IMINT  PROGRAM 


*  IMAGE  INTENSIFIED  PROGRAM 

*  OFFICE  GF  NAVAL  RESEARCH  PROGRAM 

* 


DIMENSION  Y(10)»DYU0> 

CALL  DEFINE < 1  *  6HI I DAT  ) 

READC  1  *  !  0 )  CEL*  CPC*  FREQ*A*  M*  Al< 
READ  <  1  >  1  0  >  VM*U*  UD*  S*  PASMX*  VEL 
READ <1*2 6  >  T*  GLAM  1  *  GLAM 2 
READC 1*11)  N 
WRITEC9*  12) 

12  FORM  AT  <  *EL*  ELP*  PHSE*  DPHSE  */>  • 
READC 9* 25)  EL* ELP*PHSE* DPKSE 
WRITEC  9* 1 3) 

13  FORMAT <  'IF*  Di’RT*  DT*  VDED  V) 

READ <9*  25  >  TF*  DPRT.'  DT*  VDED 
WRITE <9* 1 4) 

14  FO/-MATC  ,R*Tl*T2'/> 

READC  9*  26)  R*T1*T2 
WRITEC 9* 16) 

16  FORMAT C 'CYCMK'/) 

READC 9* 15)  CYCMX 
WRITEC 9* 1 7) 

1 7  FORMAT C ' CON VCR • / > 

READC9* 1 5)  CON VCR 

* 

10  FORMAT C  6E 1  2»  5> 

11  FORMAT < 14)  ' 

15  FORMAT CE 12* 5) 

25  FORMAT C 4E1 2. 5) 

26  FORMAT C3E 12*5) 


Figure  B2 

Page  1  of  IMINT  Program 


INITIALIZATION  OF  FLAGS  AND  SWITCHES 


"MG"  IS  INPUT  STATE  INDICATOR.  MG=1  IS  BEFORE  PULSE 
MG=  2  IS  DURING  PULSE*  MG=3  IS  AFTER  PULSE 
MG=  1 

"CYC"  IS  VOLTAGE  CYCLE  COUNTER 
CYCc 1 » 

L*=  1 
I  Gs  1 

PI=3. 14159 

"APPSW"  IS  0  FOR  FIRST  INTEGRATION  OF  FIRST  CYCLE 
APPSW»G»  * 

"TPKT"  IS  TH-E  TO  PRINT  OUTPUT  VARIABLES  . 

TPF.'T*  T 

"DTF"  IS  PERIOD  OF  DRIVING  FREO» 

DTFsTF-T 
M*  1 

"NCON**  l "  IS  SWITCH  TO  CONVERGENCE  CHECK 
"NCJNr-2"  IS  SV.HCr i  ABOUND  CONVERGENCE  CHECK  AFTER 
INTEGRATING  SXC-Ji/J  TIME  WITH  PHSE=PHSE<  SFH5E 
"NC0N«3"  Iti  S'  ITCH  COMPLETELY  AKLCViU  CONVERGENCE 
CHECKS  AMD  PHASE  ANCLE  CORRECT  TONS 
MCON= 1 

"I SAVE"  IS  FLAG  TO  SAVE  • VLL • * • T »* •TPRT * 
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FACE  . 
IMINT 


2 

PROGRAM 


570 

I  SAVE- 1 

550 

* 

. 

590 

WRITE (9, 31) 

600 

31 

FORMAT  (IX*  *CEL',9X,  'CFC  ',  9X,  'FRED  ', 

610 

WRITE <9* 30)  CEL, CPC,  FRED,  VM 

620 

WRITE <9, 32) 

630 

<M 

n 

For;  at c  1  x,  *h‘,  i  ix,  *a*,  iik,  'ai<',  iox. 

640 

WRI TEC 9, 30)  H,A,AK,V£L 

650 

WRITE <9, 33) 

660 

33 

FORMAT (IX,  'US  1  IX,  'UDS.IOX,  'S',  1  IX, 

670 

WRI  TEC  9,  30)  U,Ur.,5,R 

680 

i 

WRI TEC 9,  34) 

690 

34 

F0EMATC1X,  ‘PHSE  * , CX,  ' DPMSE ' , 7X, * TF • 

700 

WRITE C 9, 30)  PHSE,DPHSE,TF,DT 

710 

WRI TEC  9, 35) 

720 

35 

FORMAT C  IX,  * GLAM  1  ',7X,  'GLAM2',7X,  'EL 

730 

WRI TEC  9, 30 )  GLAM  1 , GLAM 2, EL, ELP 

740 

WRI TEC  9, 36) 

750 

36 

FORMAT C  1 X,  •  T1  • ,  1  OX,  ' T2  '  ) 

760 

WRI TEC 9, 30)  T1,T2 

770 

•  # 

• 

780 

30 

FORMAT C 1P4E1 2.  5) 

790 

Of  £  GA»  2*  p  I  *FR£G 

800 

DPHSEsDPHSEOPI/lSO. 

810 

PHSE*sPHSE<«FI/ 1  CO  • 

820 

■■:..r>.Cr3*0MEU  A*VH 

830 

V  E  «  7~  j  A  -  OM  C  GA/  ca.*pl>)  A 

840 

* 

' 

850 

40 

GO  TO  C45, 46), I  SAVE 

860 

45 

A1  =  T 

870 

A2=VEL 

8S0 

A3=TPRT 

890 

I  SAVE**  2 

900 

G3  TO  43 

910 

46 

TnAl 

920 

VELcA2 

930 

TPRT*»A3 

940 

♦ 

950 

43 

Pi»H*PHSE#150»/PI 

960 

* 

WRITE C 9, 840)  PPH 

970 

840  FORMAT </ 'PHASE  ANGLE  «  », FI 0.3) 

980 

WRITE C 9, 49) 

990 

49 

FORMAT C/4X,  'TIME',GX,  *VEL',6X»  'CURR 

1000 

X« 

INTEGRATION  REENTRY 

foio 

50 

V  E  L  1  c  V  E  L  5  C  0  5  C  • '  M  t:  G  A  T  +  PH  S  E  ) 

1020 

VEL4«  VM  S  7  M  C  CM  E  G  A  T + PH  SE  )  -  VEL 

1030 

GO  TO  (54,55,  56), MG 

1040 

55 

VEL6«(  1  •-EXPC-CT-T1  >/CLaM 1 > >*VEL3 

1050 

GJ  TO  58 

1060 

56 

VEL 6*  C EXP C -  C  T- 12 ) / t-L AM 2 )  >  *VEL8 

1070 

GO  TO  53 

1030 

54 

VEL6-U:-C’J:)+EL:v5?  ' 

1090 

58 

C  .•iJTI.M'Ji 

J  1  CO 

V2L«»  VEL  ,.-,V  EL  A  MATS  C  VEL  4  >  v  >R 

1  1  10 

VEL 3a- VEL «:< VEL 7  f  EXPC  AXOVEL  > 

1120 

* 

Figure  B2 
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IMINT  PROGRAM 

1130  *  BASIC  DIFFERENTIAL  EQUATION 

1140  •  VELD=(VEL1+V£L2+VEL3>/CCEL*CPC> 

1150  * 

1160  GO  TO  (60/ 1 00)/M 

1170  60  IF(T.LT. CTPRT- 1  * E- 7> >  GO  TO  80 

1180  * 

1190  CURR*(CEL*VELD>+(VEL>*(FREO**A)*H*EXP(AK*VEL> 

1200  V.'KITE(9/70>  T/VEL/  CURRi  VELD 

1210* 

1220  70  FORMAT ( 1P2E1 2*  5/ £9*  2/ El  2*5) 

1230  TPRT“TPRT+DPRT 

1240  * 

1 250  80'  CONTINUE  •  v. 

1260  GO  TO  ( 81 / 82/ 83) /MG 

1270  81  J.F C T •  LT*  ( T 1 "  1  •  E- 7 )  )  GO  TO  83 

1280  EL«ELP 

1290  MGa2 

1300  VEL8«U*(UD*EL**S> 

1310  GO  TO  50 

1320  82  IFCTvLT* (T2- 1 *E-7>  >  GO  TO  83 

1330  MG«3 

1340  '  EL«0. 

1350  VELC-U*(UD+EL**S> 

1360  GO  TO  50 

1370  83  CONTINUE 

1 3G0  * 

1390  IFCT.GE. (TF-l .E-7))  GO  TO  200 

1400  M» 2 

1410  100  CONTINUE 

1420  DY< 1 )®VELD 

1430  YC 1 ) =VEL 

1 440  * 

1450  CALL  RUNGECL/ IG/Y/DY/DT/M/T/N> 

1460  * 

1470  VELsYCl) 

1480  G3  TO  50 

1490  200  CONTINUE 

1500  GO  TO  (21 0/249/ 300  >/NC3N 

1510  v  TEST  FOR  CONVERGENCE 
1520  210  VDE«VEL 

1530  VDEAln-VDS 

1540  IF ( ABS(A3S(VDEAI >“ABS< VDED)  J-CONVCR)  300/300/230 

1550  * 

1560  *  DERIVE  NEW  PHASE  ANGLE  FOR  SECOND  CONVERGENCE  ATTEMPT 

1570  230  IK(APPSV;.E0.  !•  >  GO  TO  250 

1530  APPSV.*  1  . 

1590  PASS* 1 • 

1600  PSAVEbPHSE 

1610  PHSE"PMSE+DPH$E 

1620  N CUN * 2 

1630  GU  TO  40 

1640  240  ci:lvoe«vzl-vde 

l  650  PVDS-LELVOE/DPHSE 

1660  C’/i  TO  270 

1670  * 

1680  *  DERIVE  MEM  PHASE  ANCLE  BASED  ON  PREVIOUS  TRIES 
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PACE 

IM1NT 

1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
'  1770 
1  780- 
1790 
1800 
1810 
1820 
1630 
1840 
1850 
1860 
1870 
1680 


PROGRAM 

250  I F  < PASS • EQ • PASMX )  GO  TO  400 
PASS-PAS8+1 . 

270  COKl;aVDEAl/PVDE 
PH5E*PSAVE+C0RR 
PSAVE**PHSE 
NCONa 1 
GO  TO  40 

*  .  •  V 

300  CONTINUE 

IF<CYC.EO. CYCMX)  GO  TO  400 
IFCCYC.EO. 1 . >  NOON =3 
.  CVC---CYC-M  . 

TF=TF+DTF 

ISAVE«1 

TPRTsT 

APPSW»0* 

GO  TO  40 

* 

400  STOP 
END 


Figure  B2 

Page  4  of  IMINT  Program 
.  -  ...  ^ — ,  ,,  . 
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TABLE  BI 


Input  Parameters  for  IMINT  Program  from  Data  File 
File  Name:  IIDAT  Device  Number:  1 


Parameter 

Name 

CEL 

CPC 

FREQ 

A 

H 

AK 


VM 

U 

UD 

S 

PASMX 

VEL 

T 

GLAM1 


_ _ Comments _ ■ _ _ 

Capacitance/Cell  of  EL  [in  farads]. 

Shunt  Capacitance  around  PC  cell  [in  farads]. 

Applied  power  supply  frequency  to  the  device  [in  Hz],  . 

Exponent  of  frequency  dependence  of  EL  leakage  conductivity 
according  to:  gel  (Freq.)A. 

Dimensional  constant  of  leakage  conductivity  of  EL  [EL  leakage 
in  mhos ] . 

Exponential  coefficient  of  voltage  dependence  of  leakage 
conductivity  of  EL  according  to: 

G£L  cc  Exp (AK* Volt.),  Note:  Volt,  is  absolute  value 
of  peak  applied  voltage. 

Maximum  peak  applied  voltage  to  the  device. 

Dimensional  constant  of  PC  conductivity  [PC  conductivity  in 
mhos ] . 

Dark  conductivity  term  (U<j)  of  PC  conductivity. 

Exponent  of  steady  state  dependence  of  PC  conductivity  upon 
input  radiation. 

Upper  limit  on  number  of  searches  for  convergence  of  first 
voltage  cycle.  ,• 

Voltage  across  EL  coll  at  time  -  T. 

Time  at  which  simulation  begins  [in  seconds]. 

PC  conductivity  rise  timo  constant  used  in  form: 


GPC  1  -  Exp(-  ZlZI  )  .  (Times  in  seconds) 

GLAM1 


TABLE  BI 


Parameter 

Name 


GLAM2 


N 


Input  Parameters  for  I MINT  Program  from  Data  File  (Cont.) 
Pile  Name:  IIDAT  Device  Number:  3. 


Comments _ _ 

PC  conductivity  decay  time  constant  used  in  form: 

G__  CD  Exp(-  ) .  [Times  in  seconds] 

PC  GLAM2 

Number  of  simultaneous  first  order  differential  equations 
to  be  integrated  by  RUNGE  subroutine.  In  IMINT,  N  -  1. 


®K W IH !  PK9BIRg3|!| 


'  W  if  t*~£sV 


Parameter 

Name 


EL 


ELP 

PHSE 

DPHSE 

TF 

DPRT 

DT 

VDED 

R 

T1 

T2 

CCMX 

CONVCU 


TABLE  BIT 

Input  Parameters  for  IMINT  Program  Entered  at 
Program  Execution  Time 

_ Comments 

Input  radiation  before  input  pulse  [in  units  consistent  with 
PC  conductivity  equation]. 

Peak  value  of  square  pulse  input  radiation  [same  units  as 
EL  above]. 

Initial  guess  at  voltage  phase  angle  for  first  steady-state 
cycle  [in  degrees]. 

Increment  for  adjusting  PHSE  after  integration  of  first 
steady-state  cycle  [in  degrees]. 

Time  at  which  voltage  on  EL  repeats  its  value  at  T.  (TF-T) 
is  period  of  the  driving  frequency  [in  seconds]. 

Time  increment  between  printouts  [in  seconds]. 

Time  increment  used  per  integration  step  [in  seconds]. 

Desired  value  of  voltage  on  EL  after  integration  of  first 
steady-state  voltage  cycle  (usually  equal  to  VEL) . 

Exponent  of  input  radiation  dependence  of  PC  conductivity: 
CpQ  CO  (L)r  (Note  R  -  r-1  in  equations). 

Start  time  of  input  radiation  pulse  [in  seconds]. 

Stop  time  of  input  radiation  pulse  [in  seconds]. 

Maximum  number  of  EL  voltage  cycles  computed. 

Difference  allowed  between  EL  voltage  computed  at  end  of 
first  voltago  cycle  and  the  desired  voltago  value 
(VDED)  [in  volts]. 


a-34 


BJXJCK  DIAGRAM  OF  I  M3  NT  PROGRAM 
•  FIG  I'll  K  H3 
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Check  Status  of  Input  Radiation 
and  Adjust  Parameters  in  VEL  Derivative 
Equations  if  Necessary 


BLOCK  DIAGRAM  OF  I MINT  PROGRAM 
FIGURE  B3 
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4 


BLOCK  DIAGRAM  OF  I MI XT  PROGRAM 
FIGURE  B3 


n  -rti  i»ninmiiwHU)iin-wu»>.'. •, 


9 
_ » _ 


- - - - 

Cycle  Counter  Full?  j 

no  yes 

. _ L. .  ,  N 

I  Still  on  First  Cyclo?  I 

7  V 


STOP 


no 


yes 

1 


['NCON-3  1 


j!\^£ 


Update  Cycle  Counter 
Increment  TPRT, TF 


JL 


Save  Flag  Goes  Up 


APPSW-0 . 


BLOCK  DIAGRAM  OF  JMINT  PROGRAM 
FIGURE  B3 
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IIDAT  DATA  PILE  INPUTS  FORMAT  FOR  I  MINT  DATE: 


CEL 


CPL 


EL  Leakage  Constants 

* — : - — — -* - x 

FREQ  A  H  AK 


(6E12.5) 


S  PASMX  VEL 

_  _  _  (6E12.5) 


T  GLAM1  GLAM2 

_  _  _  (3E12.5) 


N 

(14) 


PC  Constants 

f — - — A - - 

VM  U  UD 


COMMENTS : 

TEST  DATA: 


TABLE  Bill 
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I 


i 


TTY  Inputs  for  IMINT 

EL  ELP  PHSE  DPHSE 


TF‘  DPRT  DT  VDED 


R  T1  T2 


CCMX 


(E12.5) 


CONVCR 


(E12.5) 


COMMENTS : 


TEST  DATA: 


DATE: 

(4E12.5) 

(4E12.5) 

(3E12.5) 


TABLE  BIV 
n-40 


PACE  1 


RUNGE  SUBROUTINE 

10  SUBROUTINE  f?UNC>E(L*  16*  Y*  DY*  DT*M*T»N> 

20  *  .N= NUMBER  OF  DIFFERENTIAL  EQUATIONS  (MAX.  OF  10) 

.30  *  I G= 1  INDICATES  IF  CALL  IS  FDR  FIRST  RUNGE  SLOPE 

40  *  CALC.  IN  THE  SUBROUTINE 

SO  *  T= INDEPENDENT  VARIABLE  • 

60  *  DTS INTEGRATION  INCREMENT  OF  IND.  VARIABLE  •  1 

70  *  Y= DEPENDENT  VARIABLE 

80  *  DY» DERIVATIVE  OF  Y  • 

90  *  L=2  IS  S'., ITCH  AROUND  "IG"  CHECK  AT  STATEMENT  100  (ON 

100  *  FIRST  CALL  TO  RUNGE  SLOPES)  *  ‘ 

110  *  M=2  INDICATES  TO  MAIN  PROGRAM  THAT  RUNGE  SLOPES  ARE 

120  *  ‘  BEING  COMPUTED 

130  ‘DIMENSION  DY <  1  0 ) *  Y(  1 0  )  *  F (  SO ) 

140  90  G3  TO  (100* 1 10)* L 

150  100  GO  T0( 101* 1 10)* IG 

160  101  J*1 

170  L=2 

180  DO  106  K=1*N 

190  K3=N+K 

200  106  F(K3)=Y(K> 

210  G'J  TO  90 

220  110  DO  140  K= 1  *  N 
230  K1=K 
240  K2=  J<+5vN 

250  K3-K2+N 

260  K4*K+N 

270  GO  TO  (  1  1 1  *  1  1  2*  1  1  3*  1 1  4)*  J 

280  111  F(K  1  >  =  DY(K)-:-DT 

290  Y(K)sF(i\4)  +  «5vF(Kl  ) 

300  GO  TO  140 

310  1 12  F(X2>sDY(K)*DT 
320  GO  TO  124 
330  113  F(K3)«DY(K)*DT 
340  GO  TO  134 

350  114  Y(K>  =  r(K4)  +  <F(Kl)+2.A(F(K2)+F(K3)  >+,DY(K>*DT)/6» 

360  GO  TO  140 
370  124  Y(K)s. 5*FCK2> 

380  Y(K)«Y(K)+F(K4> 

390  GO  TO  140 

400  134  Y(iO*F(X4) +F(K3) 

410  140  CONTINUE 

420  GO  T0( 1 70* 180* 170* ISO)* J 

430  170  T=T+ • 5vDi 
440  ICO  JsJ+1 

450  IF( J-4)4Q4*  404*  299 

460  299  I G= 1 

470  M® 1 

480  GO  TO  406 

490  404  I G=  2 
500  405  Lb  1 
510  406  RETURN 
520  END 

• 


Figure  B4 
RUNGE  Subroutine 
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APPENDIX  C 


UTILITY  PROGRAMS 


Cl .  Introduction 

Several  smaller  programs  were  developed  to  aid  in  the 
analysis  of  proposed  devices.  The  first  group  of  programs  are 
used  as  aids  in  extracting  useful  parameters  of  the  light-emitting 
materials  and  the  photoconductive  materials.  This  group  includes 
the  GFIT,  ELFIT,  ELDG,  and  PLDATA  programs  and  the  SLPLOT  sub¬ 
routine. 

The  GFIT  program  fits  a  non-voltage  dependent  equation  to 
resistance  vs.  input  radiation  measurements  of  photoconductor 
cells.  The  ELFIT  program  fits  the  electroluminescent  brightness 
vs.  frequency  and  voltage  measurements  to  an  equation  over  three 
consecutive  ranges  of  applied  voltage.  The  ELDG  program  creates 
pairs  of  brightness  vs.  voltage  data  from  an  equation.  The  data 
is  put  in  a  file  for  plotting  with  the  PLDATA  program.  The  PLDATA 
program  takes  data  from  a  file,  modifies  it  for  use  in  the  SLPLOT 
subroutine,  and  calls  the  subroutine.  The  SLPLOT  subroutine 
generates  a  four-cycle  semi-logarithmic  plot  of  a  maximum  of 
five  curves  of  fifty  points  each. 

The  second  group  of  utility  programs  simulate  unusual  ways 
of  driving  image  conversion  panels  and  the  stacking  or  cascading 
of  image  conversion  panels.  The  AM  and  FM  programs  use  data  files 
generated  by  the  IR  program  to  simulate  the  time  averaged  light 
output  from  image  conversion  panels  driven  by  very  simple  amplitude 
or  frequency  modulated  voltage  waveforms.  The  MASSA  program  con¬ 
denses  the  data  files  generated  by  the  IR  program  for  use  in  the 
AM  or  FM  programs.  The  CASCAD  program  uses  data  files  of  input 
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PACE  1 


RJNGE  SUBROUTINE 


10 

SUBROUTINE  RUNGE CL,  IG,  Y,  DY,  DT,M* 

T,N> 

80 

<* 

.NaNUMDEK  OF  DIFFERENTIAL  EQUATIONS  (MAX.  OF  10) 

.30 

* 

IG= 1  INDICATES  IF  CALL  IS  FOR  FIRST  RUNGS  SLOPE 

40 

* 

CALC.  IN  THE  SUBROUTINE 

SO 

* 

T* INDEPENDENT  VARIABLE 

*  ,  * 

60 

* 

DT* INTEGRATION  INC KEN ENT  BF  IND. 

VARIABLE 

70 

<« 

Y= DEPEMDEN T  VAR I ABLE 

* 

80 

* 

DY« DERIVATIVE  OF  Y  ‘ 

90 

* 

L®2  IS  SWITCH  AROUND  MIG"  CHECK  AT  STATEMENT  100  CON 

100 

* 

FIRST  CALL  TO  RUNGS  SLOPES) 

110 

* 

M«2  INDICATES  TO  MAIN  PROGRAM  THAT  RUNGE  SL0PES  ARE 

120 

* 

BEING  COMPUTED 

■ 

130 

'DIMENS  1  ON  DY  C 1 0  > , Y  C 1 0 ) , F  C  CO ) 

- 

140 

90 

GO  TO  C 100, 1 10), L 

ISO 

100  GO  TOC  101*1 10), IG 

160 

101 

J=1 

170 

L*2 

180 

DO  106  K= 1 , N 

190 

K3«N+K 

200 

106  FCK3)«Y<K> 

210 

1  10 

GO  TO  90 

220 

DO  140  K- 1 , N 

230 

K1~K 

* 

240 

}<2®}<+S*N 

250 

K3«K2+N 

• 

260 

K4~ K+N 

270 

GO  TO C  1  1  1,  1  12,  1  13,  1  14),  J 

280 

1  1  1 

FCK1 )«DY<K)*DT 

• 

290 

YCK)=FCK4)+.S*FCK1 ) 

300 

GO  TO  140 

310 

112 

FCX2)«DYCSO*DT 

320 

GO  TO  124 

330 

113 

FCK3)«DYCK)#DT 

340 

0.0  TO  134 

350 

114 

YCK)  =  FCK4)  +  CFCK1  ) +2 . *C F CK2) +F C K3 >  > +DYCK) *DT )/6. 

360 

GO  TO  140 

370 

380 

390 

1  24 

YC K) -  ♦  5 C K2) 

YCK)=YCK)+FCK4) 

GO  TO  140 

400 

134 

YCiO*FC:<4)+FCK3) 

Figure  B4 

410 

420 

140 

CONTINUE 

GO  TOC  1 70, 1  SO, 1 70, 1 80),  J 

RUNGE  Subroutine 

430 

170 

T=T+.5*DT 

440 

180 

J= J  + 1 

450 

IF  C J- 4)404,  404, 299 

460 

299  I G® 1 

470 

M=  l 

480 

GO  TO  406 

490 

404 

I  Gs2 

500 

405 

L«  1 

510 

406 

RETURN 

520 

END 
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vs.  output  o f  simulated  image  conversion  panels  to  compute  the 
light  output  of  a  cascaded  pair  of  panels.  The  MASSA2  program 
is  used  after  the  MASSA  program  to  condense  IR2  data  output  for 
use  in  CASCAD.  In  addition,  the  AM  or  FM  programs  can  be  used 
to  effectively  sort  IR  output  data  before  use  in  CASCAD. 

C2.  GFIT  Program 

This  program  fits  an  equation  of  the  form: 
g  m  conductivity  -  u(ud  +  Ts) 

where  u  and  ud  are  constants,  T  is  input  radiation,  and  S  is  a 
constant,  to  resistance  vs.  input  radiation  data  obtained  by 
measuring  individual  PC  cells  under  stimulation.  The  fit  is 
obtained  by  assuming  that  for  T  *  10“3,  g  ~  u-  udT8.  The 
natural  logarithm  of  both  sides  is  taken  and  the  resulting 
equation  is  fitted  by  means  of  the  least-squares  method  to  the 
test  data.  This  determines  S  and  U.  After  this  fit,  the  simple 
average  of  the  dark  current  data  is  used  for  the  product  u  «ud 
and  ud  is  computed. 

The  program  uses  data  from  a  stored  data  file  containing 
the  input  radiation  and  resulting  cell  resistances.  The  file  is 
called  for  at  program  execution  time.  A  zero  input  radiation 
value  signifies  the  end  of  the  data  file  and  reads  in  the  average 
(computed  manually)  dark  resistance  of  the  cells. 

Table  Cl  summarizes  the  necessary  data  inputs  and  formats 
for  the  program.  Figure  Cl  is  a  listing  of  the  GFIT  program. 
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GFIT  PROGRAM 

20  *  CONDUCTIVITY  EQUATION  fit 

40  G=U*<  UD  +  1  **S) 

60  DIMENSION  T  <  75  >  /  G<  75  >/  TL  (  75)  /  CL  <  75) 

80  DATA  S1/S2/S3/S4/  4*0./ 

100  DO  20  1*1/ 75 

120  KEADCt/2)  TCI)/ 0(1) 

140  G<I)*1»/GC1> 

160  2  FORMAT <  2E  l  6.  4) 

1 80  I E* I -  1 

200  IFC  T< I ) »LE»0. )  00  TO  30 

220  20  CONTINUE 

240  30  DO  40  J* 1/IE 

260  GL  <  J ) =  AL  0  G  <  G ( J  > ) 

280  '1 L  C  J )  •••  A  L  U  G  C  T  C  J  )  ) 

300  40  CONTINUE 

320  DO  50  KM /IE 

340  S1*S1  +GL<K)*TL<K) 

360  S2*S2+<TL(K) )**2 

380  S3*S3+GLCK) 

400  S4=  S4+  TL  <  K ) 

420  50  CONTINUE 

440  J*  N  s  1 E 

460  S*<  C(S3*S4)/FN)-S1 >/<<<S4**2)/FN)-S2> 

480  HL= <  S3- S*S4) /FN 

500  4  F'JKMATC  1P3E10.2) 

520  U=KNP<HL> 

540  UD=GCIE+1)/U 

560  WRIT E< 9/ 8) 

580  3  FUMAT(2X/  'LIGHT  IN’/SX/'G  M  £AS»S2X«'G  CALC.'/) 

600  1 1 ® I E+ 1 

620  DO  60  EM/  II 

640  GG=U*(UD  +  'i  CL)**5) 

660  '  KIT  EC  9/  4)  T < L )  /  G<L  )  /  GG 

680  60  CONTINUE 

700  WRIlE(9/5>  U/  UD/  S 

720  5  FORMA J (/ *U* •/  1  PE  10. 2// 'UD=  */ 1  PE  10. P/t  *S* ' /  1  PE  1 0. 2/ ) 

740  STOP 

760  END 
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INPUT  DATA  FORMAT  FOR  GFIT 


EQUATIONS : 

DATA  BLOCK: 
Device 

1 


G  -  U*(UD  +  T**s) 


Data  Format 


Energy  Input 


Cell  Resistance(A) 

__ _  (E16.4) 


0. 


(Maximum  of  75  lines  of  data) 


Average  Dark  Resistance 


TABLE  Cl 


C3 .  ELFIT  Program 

This  program  fits  an  equation  of  tho  form: 

B  *  Brightness  -  exp(-"^  A ), 

whore  D,  a,  and  A  are  constants;  V  is  the  applied  r.m.s.  voltage; 
and  F  is  the  applied  frequency ,  to  the  brightness  vs.  voltage  and 
frequency  data  obtained  from  sample  EL  material.  The  natural 
logarithm  of  both  sides  of  the  equation  is  taken  and  the  resulting 
equation  is  fitted  by  the  method  of  least-squares  to  the  measured 
values.  The  program  accepts  up  to  10  frequencies  and  30  voltage- 
brightness  number  pairs  per  frequency.  The  data  is  stored  in  a 
data  file  and  called  for  at  program  execution  time. 

In  practice,  it  is  best  to  first  plot  the  measurements  as 
In  B  vs  and  determine  the  Voltage  ranges  over  which  the 
resulting  curves  are  approximately  straight  lines.  Usually,  for 
thick  film  EL,  the  curves  are  straight  lines  for  all  practical 
voltages.  TFEL  curves  are  not  and  in  the  IR  program,  we  make 
provisions  for  dividing  up  the  brightness  vs.  voltage  dependence 
of  the  TFEL  into  3  voltage  ranges.  The  ELFIT  program  is  then 
used  on  the  measurements  falling  into  each  separate  range.  The 


IR  program  then  uses  the  a,  A,  and  D  parameters  only  when  the 
EL  voltage  is  in  their  range  of  validity. 

Table  cil  lists  the  data  format  for  the  stored  data  file 
used  by  ELFIT.  Figure  C2  is  a  listing  of  the  ELFIT  program. 


a-46 


PAGE  1 


ELF  IT  PROGRAM 


10 

♦ 

EL  EQUATION  FIT 

20 

* 

BaD*?F**ALRH>/EXP??A/V>**.5> 

30 

DIMENSION  CK( 10)#JJ?10)#V? 10#  30) 

40 

♦ 

# DC  10#  30 ) *  F  C 10) 

so 

DATA  FN#  S 1 #  S2#  S3#  S4#  S5#  S6#  S7#  S8/ 

80 

DO  10  I ®  1  #  10 

110 

HEAD? 1 #  2>  F ( I ) 

120 

2 

FORMAT? El 2«  3> 

130 

IF?F?I>«0.)  30#  30#  40 

140 

40 

CONTINUE 

ISO 

FI»  I 

160 

1 1  - 1 

190 

DO  20  J= 1 #  30 

200 

HEAD? 1 #  5)  V?I#J>#  B? 1#  J> 

210 

5 

FORMAT?  2E1 2*  3) 

220 

IF?V? I# J>-0. >  10# 1 0#  22 

230 

22 

CONTINUE 

240 

CF?I>sJ 

250 

Jvl?I)aJ 

260 

20 

CONTINUE 

270 

10 

CONTINUE 

280 

30 

CONTINUE 

290 

★ 

300 

DO  50  K=l# II 

310 

50 

FN=FN+CF?K) 

320 

* 

330 

♦ 

340 

DO  60  L»l# II 

350 

* 

360 

JJJ3JJ ?L ) 

370 

DO  70  MM a  1 #  JJJ 

380 

* 

390 

IF?  1 1  -  1  )  55#  55#  23 

400 

55 

S 1 SS1  +  ALUG?  Ei?L#  MM  )  > 

410 

S3*S3  +  ?V?L#MM>  >**-. 5 

420 

S7«  S7+  ?  V  ?  L  #  MM >  **- . 5 ) *  ALO  G ( B  ?  L  #  MM ) 

430 

S8sS8+?l»/V?L#MM>) 

440 

GO  TO  70 

450 

23 

CONTINUE 

460 

* 

470 

S 1  *  S 1  +  AL  0  G  ?  B  ?  L  #  MM  >  )'  / 

480 

S2aS2+AL0G?F  ?L) ) 

490 

S3aS3+?V?L#MM>>**-.5 

500 

S4aS4+?AL0G?K?L> ) >**2 

510 

S5=  S5+AL0G?  b?L#MM ) >  vALOG?  F C L  >  > 

520 

S6«S6  +  ?V?L#MM>**-.  5>*AL0G?F?L>> 

530 

S7a«7+?VCL#MM)«+-.5)*AL0G?B?L#MM> 

540 

S8»S8  +  ? l./V?L#MM>  > 

550 

70 

CONTINUE 

560 

60 

CONTINUE 

570 

* 

580 

♦ 

590 

IF? I I- 1 >  44#  44#  45 

600 

* 

610 

45 

CONTINUE 

620 

DET*FN*?-S4*S8+S6*S6>-S2*?-S2*S8+I 
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ELK II  PROGRAM 

630  DLUGe<M*C-84*SS  +  S6*i>6>-S2*<-S5*Sl3  +  S7*S6>-S3*<sr>+S6«S7 

6 AO  *  /  DET 

61i0  ALPH=<FNiK-8b*SG  +  87*S6>-Sl*<-S2*$8+S3*.S6>-S3*<SS*S7 

660  -S3*Sl»>>  /  DET 

670  Sei?A«CFN#(S^^i»7-S6*S5>-SS*(S2*S7-S3*S5>+Sl*<S8*S6- 

630  +  S3*S4>>  /  DET 

69.0  GO  TO  31 

700  * 

710  A  A  CONTINUE 
7  SO  DET«-FNvS8+S3#S3 

730  tt.UGa<»S|*S8*S3*S7>  /  DET 

740  S0PA«<FN*S7-S3*SI  )  /  DET 

700  AU’’H«0. 

760  H< 

770  if 

730  31  A  *  .*>  0  R  A  H»  .*»  0  R  A 

790  D=i£XPiDLOG> 

300  * 

810  UK' I  TEC  9>  99) 

820  99  FUNK.  AT  t  /  /  *  MOSUL TSV*  7<  •-•>> 

830  DO  DO  LL“  1  >  1 1 

840  WilU&:<9»6)  F<LL> 

bbO  6  FOKNATC//  ‘KKKO.  V,  1PK10.2/*  'VOLT.  *,5X#  ‘BRIGHTNESS*# 

860  +  SX)  'CALC*  BR*  '/> 

8 70  J J J J“ J JILL ) 

880  DO  90  Ms  1  *  JJ  J  J 

(>90  IK(II-J)  71  >  71  >  72 

900  71  LTD*"5  D/ EXP  <  <  A/V  £LL>  i'i)  )**»  5) 

910  GO  TO  73 

9 SO  72  CONTINUE 

9  3  0  » 1 0-  U  *  <  C  K  C  L  L  )  )  *  AL  r  H  )  /  E  X  /■ 1  (  (  A  /  V  (  L  L  >  M  >  )  *  *  .  5  ) 

9A0  73  CONTINUE 

9 IK)  LK1  7  EC  9*  7 )  V<LL#M>*  D<UL;M)«  Bii 

960  7  '  K'.u 6-1 A TC  1  POL  10.  R ) 

970  90  CONTINUE 

9 CO  80  CONTINUE 

990  WKIIKC  V)o>  l>/  A#  ALPH 

1000  8  KO /..'-I A  1  <//  '  D-;  '  j  l  PE  1  0.  '?./*  'A*  •• IPKI0«8/«  'ALK~  '  * E 10 •  ?//  > 

1010  STOP 

10 SO  END 
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INPUT  DATA  FORMAT  FOR  ELFIT 


EQUATION:  B  -  D*(F**ALPH)/BXP  ( (A/V)** . 5) 

DATA  BLOCK: 

Device 
1 

■ 

i 

t 

t 

i 


Data  Format 

Frequency  _ (E12.3) 

Voltage _ ,  Brightness  _ (E12.3) 


Q.  (0.  voltage  -  ready  for  new 

frequency) 

0 [* _  (0.  frequency  terminates 

data  input) 


Maximum:  10  Frequencies 

30  Volt,  Brightness/frequency 


TABLE  CII 
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C4.  SLPLOT  Subroutine 

This  subroutine  generates  a  four-cycle,  base  ten,  semi- 
logarithmic  plot  of  up  to  five  curves.  The  subroutine  uses  an 
X-array  (linear  scale),  a  Y-array  (log  scale)  and  a  variable  named 
DEC  which  determines  the  lowest  ordinate  on  the  plot.  The 
lowest  ordinate  is  10DEC. 

The  program  automatically  determines  the  end  points  of  the 
linear  scale  and  picks  a  convenient  scale  factor  for  reading 
values  from  the  plot. 

For  plotting  brightness  data  of  EL  lamps,  the  PLDATA  program 
discussed  in  the  next  section  is  used  as  a  calling  program.  Table 
CIII  summarizes  the  calling  arguments  of  the  subroutine.  Figure  C3 
is  a  listing  of  the  SLPLOT  subroutine. 


\  i 
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SUBROUTINE  SLPL0T 


20 

40 

60 

80 

♦ 

100 

120 

140 

160 

207 

180 

208 

200 

202 

220 

240 

203 

260 

280 

300 

320 

205 

340 

201 

360 

206 

380 

200 

400 

* 

420 

* 

440 

460 

4b  0 

500 

520 

540 

560 

1  1 

580 

30 

600 

31 

620 

32 

640 

20 

660 

10 

680 

♦ 

700 

* 

720 

740 

760 

780 

12 

800 

13 

820 

840 

860 

14 

880 

16 

900 

920 

940 

15 

960 

960 

1000 

61 

1020 

1040 

63 

1060 

62 

1080 

1100 

1120 

64 

208*206*206 
202*202*  201 


205*  207*  207 


SUBROUTINE  SLPL0T  <X*Y*DEC> 

DIMENSION  MCT<5)*PL0TC5>*LINEC62)>  SCALE  C  7)*XC5*50)*Y(5*5Q) 

INTEGER  ELANK*  DOT* XPLDT*  PLUS*  PLOT 

ARRAY  COMPRESSION 

DO  200  1=1*5 

MC=0 

DO  201  J= 1*50 
IF< J-C50-MC)) 

IF<X(I* J)-0.) 

MC*MC+1 
JJ=J 

X< I* JJ)  =  X< I* JJ+1  > 

Y< I* JJ>*Y<I* JJ+1 ) 

JJ=JJ+I 

IFC  <50-MC>-< JJ- 1 ) ) 

GO  TO  203 
CONTINUE 
MCT< I )=50-MC 
CONTINUE 


DEP.  VARIABLE  RANGE 
XM  IN®  X  C  I  *  1  ) 

XMAX=X< 1* I > 

DO  10  1=1*5 
MAXJ«MCT< I ) 

DO  20  J= 1  *  MAXJ 
IF<XCI* J>-0. )  20*20*11 
I K  <  X  <  I  *  J  )  -  XM I N )  30*30*31 
XMIN=X( I*  J) 

IKCXCI* J)-XMAX)  20*20*32 
XMAX=X( I*  J> 

CONTINUE 

CONTINUE 

SCALE  DEP. VARIABLE 
IX=0 

X  RAN  GE*  XM AX- XM I N 

X KSA Vs X RANGE 

I K  <  X RANGE- 60 . >  13*14*14 

X RANGE* 1 0 . ❖X RANGE 

IX=IX+1 

GO  TO  12 

IKCXKANGE-600. >  15*15*16 
X RAN GE=X RANGE/ 10* 

IX=IX-  1 
GO  TO  14 
CONTINUE 
X  D I  V=  X  RAN  GE  /  6  0  • 

IF( XDI V- 1 • >  61*61*63 
DIV*  1 » 

GO  TO  69 

IFCXDIV-2.)  62*62*64 
XDIV*2.*(10.**<-IX>> 
DIV=2. 

GO  TO  69 

IF<XDIV“5»  >  65*65*66 
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SUBROUTINE  SLPLtiT 

1140  65  XDI\'s5.*<  10.**<-IX>  ) 

1160  DIV-5. 

1180  GO  TO  69 

I  BOO  66  CONTINUE 
1 220  70  XDIVslO.*<10.**(-IX>> 

1240  l)IV~l(). 

1260  69  CONTINUE 

1 280  * 

1300  *  LOWEST  SCALE  VALUE 

1320  XS-XMIN 

1340  1 XM " 0 

1  360  80  i K ( XS- 1 • )  78*  73*  73 

1360  7  2  XS-X5U0. 

)  /him  f  ■•••-/,-  I  i 

1  420  60  TO  80 

1 440  73  JLPC XS- 10.)  75* 75* 74 

1460  74  X3-XS/10. 

1480  I  Xlvl  ~  1  XM  -  1 

1500  UU  Tb  73 

1520  }  5  XS- XS*  1  0  • 

1540  IXS--XS/S. 

1560  X5-1XS 

1380  X  S  —  a  2  *  5  • 

1600  XS-XS*<  10.**(-IXM~1  )  ) 

1620  1F(  ( XM IN- XS > • GE ♦ U  0 •  t-XDI  V >  >  X3-XS+(  10.+XDIV) 

1640  * 

1660  »•  CHOOSE  PLOT  SYMBOLS 

1600  PLANK3 1 M 

1700  DU-  1  H  -r- 

1720  XPLOTslHX 
1740  PLUS-  1H  + 

1760  PLOT ( 1 ) » 1  HI 

1  780  PLOT <  2> - t  M2 

1800  PLOTC  3)  •"  1  H3 
1820  PLOT  (4) -1114 

1840  PLOT  C  5)  -  1115 

1860  * 

1880  *  SCALE  THE  POINTS 

1900  D.J  300  1-1*3 

19  20  uAXJ-NOTt  I  ) 

19  40  L'O  301  1  *  liAXJ 

I960  XU*  J)  ~  CXO*  J)~XS)/  X  D 1  V 

I960  IP(7U*  .D.LL.O.  )  Y(1*J)"1»E“10 

2000  Y ( I  >  J ) = U.LSUlOf Y< 1* J) >*13. >-<) 8 . *DEC>  +  1 . 

2020  30)  Cv  N'UNtlE 
2040  300  CONTINUE 
2060 

2080  *  PRINT  LJKOER 

2100  SCALE ( 1 > -XS 

2120  DJ  120  1-2*  7* 1 

2140  120  SCALE < I >=5CALE< I- 1 >  +  <  1 0 .  *XL)I  V > 

2160  WR11EC9*  400) 

2180  400  KOREA  f(  ///  > 

2800  W  K I i  i .  <  9  *  121 )  SCALE 

2220  1  2 1  PU i  I  -A  TC  7  <  1  PE 8 .  :>*  2X  >  > 

2240  WHITE <9* 122) 
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SUBROUTINE  SLPL0T 

2260  122  F0RMATC  7C  •  1  '  /  9X)  ) 

2280  D0  123  1=1/61 

2300  123  LINE< I )=PLUS 
2320  DO  124  1-1/61/10 

2340  124  L INE  C  I )  =XPLOT 

2360  WRITE (9/ 125) (LIME C I DUM >/ IDUM=1/61 ) 

2380  1 2b  F0  KM AT C  6 1  A 1  ) 

2400  *  BLANK  THE  LINE 

2420  DO  126  1=1/61 

2440  126  LI NEC  I )  =  BLANK 
2460  * 

2480  *  PLOT  THE  POINTS 

2500  IDEC=DEC+4. 

2520  ID0T=O 

2540  IXPL0T*O 

2560  L INEC  62 )  =  0 • 

2580  DO  310  N= 1/71 

2600  LINEC1)=PLUS 

2620  I DOT- ID0T+ 1 

2640  IFC IDWT.NE.6)  GO  TO  315 

2660  ID0T=0 

2680  LINE<l)=DOT 

2700  IXPL0T=  I  XPLUT+ l 

2720  IF C IXPLOT.NE. 3)  GO  TO  315 

2740  IXPLUT=0 

2760  LINE( 1 )=XPLUT 

2780  LINEC61 >=XPLUT 

2800  IDEC= I  DEC* 1 

2820  L INEC  62 )  =  I  DEC 

2840  315  IPY=73-N 

2860  DO  311  1  =  1/5 

2880  MAXJ-MCTCI) 

2900  DO  312  J= 1 / MAX J 

2920  IY=YCI/ J)*l .001 

2940  IX=XCI/J) 

2960  IKCIY-73)  316/317/317 

2980  31 7  YC 1/ J)=0. 

3000  I Y=0 

3020  316  IFCIY-IPY)  312/313/313 
3040  313  L1NECIX+1 )=PL0TCI > 

3060  YCI/J)=0, 

3080  312  CONTINUE  - 

3100  311  CONTINUE 

3120  1 0=6 1  FIGURE  C3 

3140  DO  314  I DUM= 1/61 

3160  IF  CL1NEC  IO.WE# PLANK)  GO  TO  384  - —  ■ 

3180  314  I  C= I C-  1 
3200  384  I C= I C+ 1 

3220  WKXTEC9/380)  CL  INK  CN DUM )/ NDUM* 1 / I C> 

3240  380  FORMAT  C61A1 / * l »0E  */ 12) 

3260  DO  382  11=1/61 

3280  382  LINEC I I )= BLANK 
3300  310  CONTINUE 
3320  DU  383  1=1/61 

3340  383  L  INF C 1 >=PLUS 
3360  DU  385  1=1/61/10 
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SUBROUTINE  SLPLOT 


3380 

3400 

3420 

3440 

3460 

3480 

3500 


385  LINE< I >=XPLOT 

WHITE ( 9* 1 25) CLINE ( I  DUN >  ,  I  DUM= 
WRITE ( 9# 1 22) 

WKITEC9* 1 21 >  SCALE 
WRITE <9* 400) 

RETURN 

END 


*61  ) 


DATA  BLOCK  FOR  SLPLOT  SUBROUTINE 

<X,  Y,  DEC) 

Dimensions  of  Arrays: 

X( 5, 50) 

Y(5, 50) 

X  direction  is  linear  scale 
Y  direction  is  Log1Q  scale 
Maximum  of  5  curves  plotted 
DEC:  Bottom  ordinate  of  the  Plot  is  10DEC 


TABLE  CIIX 
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C5.  PLDATA  Program 

PLDATA  converts  raw  measurements  of  brightness  of  EL  vs. 
voltage  at  a  fixed  frequency  into  a  data  block  for  plotting  by 
the  SLPLOT  subroutine.  The  program  lets  you  choose  at  execution 
time  whether  to  plot  vs.  voltage  or  the  inverse  square  root  of 
the  voltage.  This  second  form  allows  the  operator  to  determine 
ranges  of  applied  voltage  over  which  the  ELFIT  program  can 
provide  accurate  results. 

Table  CIV  summarizes  the  format  of  the  data  file  used  with 
the  program.  Figure  C4  is  a  listing  of  the  PLDATA  program.  It 
must  be  loaded  with  the  SLPLOT  subroutine. 
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HA  GE 


1 


HLDATA  PROGRAM 


20 

* 

PLOT  DATA  WITH  SLPL3T 

40 

DIMENSION  XC5# 50># Y<5# 50) 

60 

DATA  X/250*0./# Y/250+0./ 

80 

DO  10  I«l#5 

100 

DJ  20  J= 1 #  50 

120 

READC 1  #  1  )  X<I# J># Y< I# J> 

140 

IF<X( I# J) »LE« 0» )  GO  TO  21 

160 

20 

CONTINUE 

180 

21 

CONTINUE 

200 

10 

CONTINUE 

220 

READC 1  #  g>  DEC 

240 

1 

F0KMATC2E1 6.  3) 

260 

£ 

FORMAT <E  1  6» £> 

280 

WRITE <9#  :o 

300 

3 

FORMAT </// ’STYLE  OF  PLOT'/ 

320 

+ 

'1=  X- VALUES  AS  READ  V 

340 

+ 

'2-  X=1/S0RTCX>V 

360 

+ 

) 

380 

KEADC  9#  4>  IS 

400 

4 

FOKMATCIl) 

420 

GO  TO  <30#  40  ># IS 

440 

40 

DO  50  I -  1  >  5 

460 

DJ  60  1 #  50 

480 

IFC  X< I# J> .NE*  0. >  X< I#  J>«  1 • /SORTC  XC  I 

500 

60 

CONTINUE 

520 

50 

CONTINUE 

540 

GO  TO  100 

560 

30  CONTINUE 

580 

100 

CALL  SLPLO’IC  X#  Y#  DEC  > 

600 

IF< I EOF  Cl))  102# 102# 103 

620 

102 

KEADC 1, 101 > 

640  101  FORMAT ( * 

660  WRITE <9> 10! ) 

600  WRITE <  9# 200) 

700  200  FORMAT C /// ) 
720  103  S’UP 
740  END 


FIGURE  C4 
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INPUT  DATA 


Device 

1 


FORMAT  FOR  PLDATA  (Uses  SLPLOT  Subx'outine) 


DATA  FORMAT: 

Linear  Scale  Log  Scale 

X _ ,  Y _ (2  E16.3) 


i  » 


<  • 

i  • 

0«  (terminates  a  curve) 

0. 

0. 

0, 

_ 0.  (After  five  zero.  X  inputs, 

ready  for  DEC) 

DEC  _  (E16. 2)  (Lowest  ordinate  plotted  -  10uliU) 


(Optional  30  Character  Line  for  Data  Identification) 


Beware  of  Blank  Records  in  File!! 


TABLE  CIV 
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C6.  ELPG  Program 

This  program  generates  a  set  of  EL  brightness  vs.  voltage 
and  frequency  values  computed  from  the  equation: 

Brightness  -  D»Fa  exp(-V 7  ) 

where  D,  a,  and  A  are  constants,  F  is  the  applied  frequency,  and 
V  is  the  applied  voltage.  This  program  is  useful  in  checking  the 
results  of  the  ELFIT  program.  The  resultant  data  file  is  ready 
for  use  by  the  PLDATA  program. 

The  program  asks  for  the  D,  a  and  A  values  along  with  their 
ranges  of  validity  at  execution  time.  All  data  is  entered  at 
the  keyboard.  The  program  will  accept  five  voltage  ranges  but 
the  IR  program  can  only  use  three  ranges. 

Table  CV  lists  the  formats  for  the  input  data.  Figure  C5 
is  a  listing  of  the  ELDG  program. 


FACE  I 


ELDG  PROGRAM 


20 

* 

EL  DATA  GENERATOR  FOR  PLDATA 

40 

DIMENSION  F? 10>*A?5>* VL ? 5 > * VU? 5>* DV? 5 

60 

DO  10  1=1*5 

80 

WRITE? 9*  l  > 

100 

1 

FORMAT? ’VOLTAGE  RANGE)  VL#VU’/> 

120 

READ?  9*  2)  VL  ? I >  *  VU ( I ) 

140 

2 

FORMAT? HE  I 0.2) 

160 

IF?  VL?  I  )-Q»  )  30*  30*  20 

180 

20 

1 1  =  J 

200 

WRITE? 9*3) 

220 

3 

FORMAT? 'DV  *  AV  *  ALPHA'/) 

240 

READ? 9*  4)  DV? I ) * AV? I ) * A? I > 

260 

4 

FORMAT? 3E1 0.2) 

280 

10 

CONTINUE 

300 

30 

continue 

310 

2ER0=O. 

320 

* 

340 

WRITE? 9*  5) 

360 

5 

FORMAT?  ’FREQUENCIES  ' /  ) 

380 

DO  40  J=1 *  10 

400 

READ?  9*  6 )  F?  J ) 

420 

6 

FORMAT? El  0.2) 

440 

IF?F?J)-0.)  50*50*60 

460 

60 

JJ=J 

480 

40 

CONTINUE 

SOO 

50 

CONTINUE 

blO 

WRITE? 9* 7) 

51  1 

7 

FORMAT?/ ’EXPONENT  OF  SMALLEST  ORDINATE 

512 

READ? 9*  to  DEC 

513 

8 

FORMAT?  F 1 0* 1 ) 

520 

♦ 

S30  DO  80  L=1*JJ 

I>40  DO  70  K=l*  II 

660  C= 1 • 

680  V=VL?K> 

700  INTV5=VL?K>/5. 

720  VS*lNTVb 

740  VM5*V5*b. 

760  90  CONTINUE 

B40  i F?  V“ VU  ? X  > )  21*22*22 

860  2 2  V=VU?M> 

800  0  200. 

900  21  CONTINUE 

920  ELUGa  S’.'RT?  AV?  K)/V> 

9  AO  FLD(>=A?K)*ALOG?F?L> > 

960  DVLOG=ALD(i?DV?K>  > 

980  bLOG*>bVLOG+KLGG- RLOG 

1000  B* EXP?  FLOG) 

1020  WRITE (1*9)  V*  8 

1040  9  F 0 KM  A 1  ?  HE  1 6  •  3  ) 

1080  IK (C.E,’.  I .  >  VsVMb  +  5. 

1100  IF<  C*  GT» 1 • )  V  *  V  +  5 » 

1120  C=CM. 

1140  IF?C.GE. 200. )  GO  TO  70 

1160  GO  TO  90 

1180  70  CONTINUE  _ 

CxCUKE  C3  [ 
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ELDG  P HUG HAM 


1188 

1190 

1192 

1193 

1194 
1195 

1196 

1197 
1  200 
1220 


WRITE? 1, 9 >  2ER0#ZEK0 
80  CONTINUE 

IF?JJ»£G»S>  GW  TW  99 
DO  98  KK=JJ,4,  1 
WRITE? 1#  9)  2KK0* ZEK0 
98  CONTINUE 
99  WRITE? 1/ 11)  DEC 
11  FORMAT?  E 1 6  •  2) 

STOP 

END 


FIGURE  C5 
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INPUT  DATA  FORMATS  FOR  ELDG 


EQUATION:  B  -  (F(L)  **A(K) ) *DV(K)/SQRT(AV(K)/V) 


DATA  BLOCK: 


Device 


TTY 

(as  asked  for) 


£aaaftt,B 

VL 

,  vu 

(E10.2) 

DV 

,  AV 

,  ALPHA 

VL 

-  .  vu 

DV 

.  ,  AV 

,  ALPHA 

(E10.2) 


_ 0. 

Frequencies 


(0  Voltage  terminates 
inputs) 

(E10.2) 


parameter 


- (0  frequency  terminates  data  input) 

Exponent  of  (El 6. 2) 

Smallest  ~  ’ 

Ordinate 


Maximum  of  5  voltage  ranges 
Maximum  of  10  frequencies 

VL  <  Voltage  Range  <  VU 


TABLE  CV 


C7.  FMj,  AM  and  MASSA  Programs 

The  MASSA  program  condenses  the  data  file  created  by  the  IR 
program  into  a  more  compact  form  for  use  with  the  AM  and  FM  programs. 
Figure  CG  lists  the  MASSA  program  .  It  calls  for  the  input  data 
file  at  program  execution  time. 


The  FM  program  uses  the  brightness  vs.  radiation  input  data 
generated  by  the  IR  program  and  condensed  by  the  MASSA  program  to 
approximate  the  operation  of  an  image  conversion  panel  driven  by 
several  applied  frequencies  in  rotation.  The  simulation  is  that 
a  panel  is  driven  at  a  constant  voltage  and  one  frequency  (F^)  for 
a  short  time  (t^),  then  driven  at  another  frequency  (F2)  for  a 
short  time  (t2),  and  so  on  up  to  Fn  for  time  tn  and  then  driven 


again  at  F-i ,  etc.  If  the  sum  I!)  t*  is  less  than  the  time  for  the 

i-1 

eye  to  see  flicker  in  the  brightness  of  the  panel,  the  viewer 

perceives  a  simple  average  brightness  of  the  panel.  The  FM 

program  sorts  through  the  response  curves  generated  by  the  IR 

program  and  weights  the  output  brightnesses  for  a  selected  frequency 

n 

by  the  fraction  tj/£  t^  for  that  frequency.  The  program  will  accept 

Ji-1 

up  to  20  desired  frequencies  but  seldom  is  it  worthwhile  to  construct 
a  data  base  large  enough  for  use  with  that  many  frequencies. 

Figure  C7  is  a  listing  of  the  FM  program.  Table  CVI  is  a 
summary  of  the  data  file  format  for  the  FM  program.  The  FM  program 
may  also  be  used  to  effectively  sort  and  condense  data  for  use  in 
the  CASCAD  program.  The  FM  program  must  be  loaded  with  the  PLOT 
subroutine. 

The  AM  program  works  in  an  Identical  manner  to  the  FM  program 
only  it  computes  a  weighted  average  of  different  voltages  applied 
for  short  time  periods  at  one  frequency.  Figure  C8  is  a  listing 
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Of  the  program.  Table  OVXI describes  the  formate  of  the  data  m6 
This  program  may  also  bo  used  to  sort  data  for  CASCAD.  The  AM 
program  must  be  loaded  with  the  PLOT  subroutine. 
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1 

MAS3A 

PKdGKAM 

20 

* 

massage  light  amp 

DATA 

30 

IFLAG= 1 

32 

KEADC  1*4)  LP 

34 

WHI TE  <  2*  4 )  LP 

36 

'  4 

FUHMATC I2> 

40 

10 

HEAD  <  1  *  1  )  B*  T*  V*  F 

60 

1 

FURMAK  4E  1  2.4) 

80 

12 

1  =  0 

100 

BS=  Q 

1  20 

TS*T 

140 

vs=v 

160 

fc'S=K 

180 

IvKI  1  E  (  2*  2)  V*  F 

200 

2 

FORMAT ( 2c 1 2.4) 

210 

* 

220 

20 

KEAD<1*1>  B»  T*  V*  F 

240 

IF  C  I EOF < 1 ) )  30*  30 * 

100 

260 

30 

1KF-K6)  40*80*  40 

230 

40 

IFCI-19)  60*70*70 

300 

60 

DJ  80  J=l* 1 8* 1 

320 

l\  K I  Tfc:  <  2*  3 )  BS 

340 

3 

FUkMAT<E12.  4) 

360 

80 

CUN flMUE 

370 

GO  10  <70*110)*  I FLAG 

380 

70 

GO  It)  12 

400 

* 

420 

SO 

'.-.KITE <2*  3)  13S 

440 

1  =  1  +  1 

460 

bs=b 

4b0 

1  8=T 

800 

Vb=  V 

b20 

F8=F 

b40 

GO  10  20 

860 

♦ 

880  100  IFLAG=2 
8*0  CiJ  U  40 

8*8  110  31  UP  fc.UK 
600  END 


FIGURE  C6 
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FM  PROGRAM 

20  *  FM  DRIVEN  LIGHT  AMP 

40  * 

60  DIMENSION  FDC  80) >FKA<  20) / B0<  19/5)*VS(5)*XX<  19*5) 

GO  COMMON  XX/  BO*  LP 

100  DATA  VS/5*0./ 

120  Dd  23  1*1*19 

140  Du  24  J=  1/5 

160  XX(1/J>  =  0. 

180  Bi)CI*J>»0* 

200  24  CONTINUE 

220  23  CONTINUE 

240  tvRI  Tr.'f  9/4) 

260  A  FORMAT </  'CREATE  DATA  FILE?  0=N0/  1  =  YESV> 

280  READ  (9/5)  I  St;  IT 

300  5  FORMAT ( I  1 ) 

320  1 SW  IIs  I 3W1T+ 1 

340  WRIT  EC  9/ 99) 

360  99  FORMAT  </  Hi  ANT  PLOT?  0=N0/  1  «YES  V  > 

360  READ (9/ 5)  I  SUITS 

400  i  St.'  1 12=  I  SW I  T2<- 1 

420  Will  TEC  9*  2> 

440  2  F  U  KM  A  T  <  /  1 F  HE  0  •  /FRACTIONS  V  ’TERMINATE  WITH  O.V> 

460  I GUUD= I 

480  DO  10  1=1/20 

500  READ <9* 1 )  FDC I ) / FKA< 1 ) 

520  1  FORMAT C8E12# 4) 

540  1FCFD(I)>  11/12/11 

560  1  1  IC--I 

560  10  CONTINUE 

600  12  CONTINUE 

620  FS=FD(1) 

640  * 

660  READ C 1 / 7 )  LP 

680  7  FORMAT C 12) 

700  IFCJSVsIT.F.0.2)  KRlTE<2/7>  LP 

720  20  READC1/1)  V'/F 

740  DO  30  1  =  1/ IC 

760  I F ( A BG ( FD ( I > “F  )  -  1  •  )  21/21/30 

760  21  1  l;'J  JD  =  8 

800  1 SAv E“ I 

820  GO  TO  22 

640  30  CONTINUE 

860  I 6JUU  = 1 

880  22  CONTINUE 

900  GO  TO  (40/60)/  I  GOOD 

920  40  DO  -60  0=1/19 

940  .  KEAD(  1/3) 

960  3  FORMAT <  E 1  2»  4) 

980  50  CONTINUE 

1000  IFCIEOFU))  20/80/80 

1020  * 

1040  60  DO  69  K= 1 / 5 

1060  IF<A»;S<V-VSCK>)-1  .  >  66/66/69 

1080  69  CONTINUE 

1100  DO  68  K= 1 / 5 
1120  IFCVSOO)  68/67/68 


FIGURE  C7 
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FM  PROGRAM 


1  1  40 

6S 

CONTINUE 

1160 

1 GOODa  1 

1  I  80 

GO  TO  22 

1200 

67 

VSCK)=V 

1220 

66 

CONTINUE 

1240 

DO  70  1=1/19 

1260 

KEA0< 1/3)  B 

1260 

B0 (  I  /  K  >  =  00 ( I / K )  +  ( FRA ( I SAVE )  *B ) 

1300 

70 

CONTINUE 

1320 

F3=K 

1340 

IFUEOF(l))  20/20/80 

1360 

80 

CONTINUE 

1380 

GO  10(91/92)/  I  CWn 

1400 

92 

CONTINUE 

1420 

DO  93  K=  1  /  5 

1440 

DO  94  J=l  /  19 

1460 

IF<f30(J/K>)  93/93/95 

1  480 

95 

ERI  '1E<2/  6)  BJ(J/K) 

1500 

6 

FORMAT (El  2.4) 

1520 

94 

continue 

1  540 

93 

CONI 1NUE 

1  560 

91 

GO  JO  <  9  7/  96)/  ISV.I  T2 

1560 

96 

CALL  '  OPL 

1600 

97 

CONTINUE 

1620 

STOP 

1640 

END 

FIGURE  C7 
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INPUT  DATA  FORMATS 

Device ; 

TTY 

(as  asked  for) 


1 


FOR  FM 

Data  Format: 

Frequency  ,  Fraction  (E12.4) 


_  0 «  (maximum  of  20) 

Voltage _ ,  Frequency  (E12.4) 

Brightness _  (E12.4) 


19  values 


Voltage _ ,  Frequency 

Brightness 


etc . 


TABLE  CVI 
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AM  PROGRAM 


20 

40 

60 

80 

100 

120 

140 

160 

180 

200 

220 

240 

260 

280 

300 

320 

340 

360 

380 

400 

420 

440 

460 


*  AM  DRIVEN  LIGHT  AMP 

* 

DATA  FS/5*0./ 

DO  23  1  =  |  ,  1  9 
DO  24  Js i , 5 
XX*I, J)-0. 

80*1, J>=0. 

24  CONTINUE 
23  CONTINUE 
WRITE <9,  4) 

4  F  0  KM  AT  <  /  *  CHEAT  E  DATA  FILE?  0=N0,  i  =  YES*/) 

READ*  9,5)  ISWIT 
ISWliBl5lvIT+l 
5  FORMAT < 1 1 ) 

WHITE < 9x 99 > 

99  FORMAT  (/  *  WANT  PLOT?  0=NO,  lsYESV) 

READ <9, 5)  I SW IT 2 
ISWM  2~  I  OWI  i  2+1 
W  k  1  I E  t  9  x  2  ) 

2  FORMAT  < / ' VOLTAGES*  FRACTIONS V *  TERMINATE  KITH  0. »/) 

1  GOODS  1  r  X  *n  u.  f  ! 


480 

K  =  0 

500 

DO  10  1=1x20 

520 

READ*  9,1)  VD(  1  )  x FRA (  I  ) 

540 

1 

FORMAT* 2EI 2.4) 

560 

IFCVD(I))  11x12,11 

580 

1  1 

IC=I 

600 

10 

CONTINUE 

620 

12 

CONTINUE 

640 

Vii=VD*  1  ) 

660 

* 

680 

READ* 1, 7)  LP 

700 

7 

FORMAT* 12) 

720 

IF*  1:>V;I  I*E0*2)  WRITE*2x7)  LP 

740 

20 

READ* lxl)  V,F 

760 

DO  30  I  =  1 ,  I  C 

780 

IF*  APR*  VI)  <  I  >-  V  )  -  1  .  )  21,21x30 

800 

21 

I  GOOD -  2 

820 

I SAVEs I 

840 

80  JO  22 

860 

30 

CUNT INUE 

880 

lOOivDs  1 

900 

22 

CONTINUE 

920 

80  TO  *40x60),  I  GOO D 

940 

40 

DO  50  Jsl, 19 

960 

READ* 1,3) 

980 

3 

1*  DMATCEt  2*  4) 

1000 

50 

CON  T INUE 

1020 

IFdEOFtl))  20*  20,  60 

1040 

* 

1060 

60 

DO  69  K= 1 , 5 

1080 

11*  *ADO<l*“f  S(K))-1  *)  66x66,69 

1100 

69 

CONI INUE 

1120 

DO  68  K=|*5 

FIGURE  C8 
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AM  PKOGHAM 

1140  IF<F S<K>)  6B/67/68 

1160  68  CONTINUE 

1180  IG0UI)=1 

1200  GO  10  22 

1220  67  FS(K  >=F 

1240  66  CONTINUE 

1260  UJ  70  1=1/19 

1280  KEALX1/3)  8 

1  300  DO  (  1  *  K  )  =  BO  <  I  *  K  )  *  <  F  i<A  <  ISAVE>*B> 

1320  70  CONTINUE 

1340  VS=V 

1360  1F< IE0F< 1 ) >  20/20/80 

1300  i j 0  CON  1 INUE 
1 400  (  J  TO  <91/ 92) / I  SO 1 T 

1420  92  CONTINUE 

1440  DU  93  K=  1  /  Sj 

1 460  DO  94  J= 1 / 1 9 

1480  I F  ( B J ( J/ K >  >  93/  93/  95 

1500  95  OKI TE< 2/ 6)  BO<J/K> 

1  520  6  FORMAT < E 1 «• 4) 

1540  94  CONTINUE 

1560  93  CONTINUE 

1  580  9 1  GO  TO  <97/  967/  I  S'..  I  T2 

1600  96  CALL  LOCPL 

1620  97  CONTINUE 

1640  STOP 

1660  END 


FIGURE  C8 
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INPUT  DATA  FORMATS 

Device : 

TTY 

(as  asked  for) 


FOR  AM 


Data  Format: 

Voltage  ,  Fraction _ (E12.4) 


_ 0 .  (maximum  of  20) 

Voltage _ Frequency _ (E12.4) 

Brightness _ (E12.4) 


19  values 

t 

Voltage _ ,  Frequency 

Brightness _ ___ 


etc . 


TABLE  CVII 
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C8 .  CASCAD  Program  and  MASSA  2  Program 


The  MASSA  2  program  converts  the  data  files  generated  by  IR 
and  MASSA  into  a  form  compatible  with  the  CASCAD  program.  The  IR 
program  can  generate  up  to  25  output  brightness  curves  for  different 
voltage  and  frequency  combinations.  The  CASCAD  program  will  only 
plot  data  for  5  possible  output  combinations  of  the  second  stage 
panel.  The  operator  should  edit  his  data  files  to  make  sure  that 
the  files  used  in  device  2  contain  no  more  than  5  curves.  The  AM 
or  FM  program  can  be  used  to  sort  the  files  generated  by  IR  and 
MASSA.  For  instance,  if  the  operator  using  AM  selects  100%  of 
one  voltage  to  be  used,  the  output  data  file  will  contain  brightness 
curves  at  only  that  operating  voltage  and  all  frequencies  that  were 
available. 

The  CASCAD  program  takes  a  set  of  19  output  brightnesses  from 
the  first  image  panel  and  uses  these  brightnesses  with  linear 
interpolations  from  all  sets  of  output  brightnesses  for  the  second 
state  image  panel  to  generate  an  approximate  simulation  for  a  two- 
stage  image  converter.  The  program  permutes  all  curves  of  the 
second  stage  panel  on  each  curve  of  the  first  stage  panel. 

Either  or  both  of  the  simulated  image  conversion  panels  may 
have  been  AM  or  FM  simulated  panels.  At  this  stage  in  the  simulations, 
operator  discretion  is  mandatory  or  one  may  drown  in  output  curves. 

Figure  C9  is  a  listing  of  the  MASSA  2  program.  Figure  CIO 
is  a  listing  of  the  CASCAD  program.  Tablccvni  is  a  summary  of  input 
data  formats  for  CASCAD. 
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MASSA2  PMGRAM 
20  *  MASSA2 

AO  *  PREPARES  .MASSAGED  DATA  FOR  CASCAD 
60  DIMENSION  B < 19) 

80  KEALK  1,1)  LP 

100  WK1 1E< 2, 1 >  LP 

120  1  FORMAT? 12) 

1  AO  2  FORMAT  <212.4) 

160  3  FORMAT? 2E 1 R.  A) 

180  30  READ? 1 , 3)V,F 

200  DO  1 0  I  *  1 »  1 y 

220  READ?  1,2)  B?  I  > 

2 AO  10  CONTINUE 
260  Uv)  HO  I=|,I9 

280  l\‘KI  TK?  2*  2)  fi<  I  > 

300  20  CONTINUE 

320  IT  <  U.JK?  1 ) >  30, 30, AO 

3 AO  AO  i.T'Jp  EOF 

360  END 


FIGURE  C9 
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PAUL 


1 


CAbCAD  PROGRAM 


20 

♦ 

CAbCAD 

40 

* 

CCMPUTfclS  OUTPUT  OF  CASCADED  LIGHT  AMPS 

00 

DIMENSION  H 1  d  9  )*  02 C  1  9  >  *  T  €  1  9  )  *  ID <  1  9*  b > *  XX C  1  9*  5 > 

AO 

COMMON  XX*  13J*  LP 

100 

DATA  1  <  1  >*  T<2>*  TC3)/1  .fc>6*  2.  E- 6*  b .  E- 6/ 

1  9.0 

DATA  T  d9)/l  ./*  J/0/ 

1  AO 

DU  12  1=1*19 

160 

DU  13  J  =  1  *  S 

1  AO 

XX  ( 1  *  J )  =  0 • 

800 

13 

CUM 7 INUE 

9.9.0 

12 

CONTINUE 

240 

DJ  10  1=4*16*3 

860 

I  d  )  =  1  d-3)*10. 

260 

1  C  1  +  I  )  *  1  <  I  -  8  )  * 1  0 . 

300 

T(I+2>«1<I-I  )  M0. 

380 

10 

CUM  UNDE 

340 

* 

360 

WMTL<9*  1  > 

360 

1 

FORMAT </ 'CREATE  DATA  FILE?  0*N0*  l = YES  */  ) 

400 

KLAD<y*2>  I  SUIT 

A  9.0 

2 

F  U/.MAl  C  I  1  ) 

440 

i  sv.  n  =  i  su  it  +  i 

460 

WRITE  <9*  3) 

400 

3 

fuRMAI  </ 'BEGINNING  INPUT  RADIATION  V) 

bOO 

READ (1*6)  LP 

b20 

WF! IT  E (  9 *  7 )  LP 

b40 

REAUC  2*  6 )  LP 

b60 

WRIT E< 9* 7)  LP 

bOO 

IF<  IbKI  I  .LO.  2)  V.'Rl TL(  3*  6>  LP 

600 

6 

FORMAT (12) 

6  9.0 

7 

FORMAT  d OX*  '  1 •£  *  12) 

6  40 

4 

h OKMAT  (El 2. 4> 

6  60 

Tn«10.**LP 

680 

+ 

700 

DJ  20  1=1*19 

7  9.0 

T  d  )  =  1  (  I  >  *  (  T;3/ 1  •  E-  6  ) 

7  40 

20 

CUnT  I WD  1C 

760 

* 

7  80 

30 

IF<  IfDp<  1  >.t;o.  1  >  GU  to  101 

800 

Du  31  1  =  1*1 9 

8  80 

RKAD(  1  *  b)  hl(I) 

840 

31 

CU.v'l  INUE 

860 

b 

FORMA  1 (812.4) 

880 

3b 

CUM  1 I WUR 

900 

DU  32  1  =  1*19 

‘>80 

KKADC  2*  b)  F3;’d> 

9  AO 

32 

CUNT  ImUL 

960 

J=J*l 

960 

N=2 

1000 

40 

DU  bO  1=1*19 

1 0  80 

Du  60  K=tM*  I  9*  1 

1040 

IK<MICI>-1(K>>  70*70*80 

1060 

80 

IF  (  K“  19)  60*90*90 

1060 

60 

CDNJ  1NUL 

1100 

90 

DJ( 1  * J)«02(K> 

1  l  BO 

CJ  ID  bO 

FIGURE  CIO 

a-7'1 


HAUL* 

CASCAD 

11  AO 
1160 
1180 
1200 
1220 
12A0 
1  260 
1280 
1 300 
1320 
I  3  AO 
1  360 
1380 
1  400 
1  A  20 
1  A  40 
1  A60 
1  A80 
1500 
1520 
1  5A0 
1  560 


PM  GRAM  ' 

70  NsK 

B0<I>  J>»D2<K-1  >  +  <Cn2CK)-B2(K-  l  >>*CB1  (I)-T(K-l  >  > 
+  /<TCK>-TCK- 1 )  >  > 

IF  C  BO ( 1 i J ) «Lt • 0 • )  B3< l» J)a 1 . E- 10 
50  CkMUNUC 

IF  < ItOF ( 2) • t0» 1 )  0J  10  100 
00  TO  35 

* 

* 

100  GO  10  <  103#  102)  *  151.  IT 

102  F 81 12(3* A)  (<B0(L#M)/L*1/ 19>#i*l*l> J» l > 

103  CALL  L0GFL 

DO  25  1  =  1>.J*1 
DO  26  K*  1  *  1  y 
B0<  K>  I  >  =  0. 

26  CUN  TIiMUL 
25  CONTINUE 
J=0 

REWIND  2 
G td  TO  30 

101  STOP 
END 


FIGURE  CIO 


INPUT  DATA  FORMAT  FOR  CASCAD 


Device: 

1 

(First 


2 

(Second 


Data  Format: 

LP _ (1 2)  (Beginning  input  Radiation 

Stage)  10^) 

Output  Brightness _ (E12.4) 


t 


19  values 

i 


Brightnesses  are  in  groups  of  19 

Stage)  LP _  (12  ) 

Output  Brightness  _ 


19  values 


(Same  format  as  for  device  1) 

(Maximum  of  5  groups  of  Brightnesses) 


TABLE  CVII I 


